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Abstract 

The nucleoli structure functions probed in deep-inelastic scattering at large virtualities form an 
important tool to test Quantum Chromdynamics (QCD) through precision measurements of the 
strong coupling constant a s (M^) and the different parton distribution functions. The exact knowl- 
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edge of these quantities is also of importance for all precision measurements at hadron colliders. 



D . During the last two decades very significant progress has been made in performing precision calcu- 

lations. We review the theoretical status reached for both unpolarized and polarized lepton-hadron 
scattering based on perturbative QCD. 

>: 

22" 1 Introduction 
o 

. Matter consists of regular structures at microscopic distances, which exhibit themselves at the crys- 
O! talline, molecular, and atomic levels pp. The discovery of a,/3 and 7 radioactivity [2] provided new 
^ ■ natural probes beyond the visible spectrum of light and X-rays to resolve even smaller structures of 
■ • ■ matter. In 1911 E. Rutherford discovered the atomic nucleus of a size much smaller than that of atoms 
• rH ■ through scattering of a-particles at gold [3]. Herewith the picture of matter at small distances changed 
^ ! dramatically rising the question for further sub-structures. The composite nature of nuclei could be 
03 I explained after Chadwick's jl] discovery of the neutron and Yukawa's model for nuclear forces [5]. An- 
other important discovery was made by Frisch and Stern in 1933 measuring the anomalous magnetic 
moment of the proton with a different value from that of point-particles, like electrons [6]. Later in 
1939 Alvarez and Bloch measured the anomalous magnetic moment of the neutron [7], both of which 
constituted first evidence on the compositeness of nucleons. The current values of the nucleon magnetic 
moments are [5] 

Hp = 2.792847356 ± 0.000000023 fi N , /i n = -1.9130427 ± 0.0000005 fx N , (1.1) 

with fiN = eh/2m p the nuclear magneton. 

During the 1950ies the Hofstadter experiments [9] operated at virtualities being large enough to re- 
veal the charge distribution inside nucleons, which is illustrated in Figure HJ A positive core distribution 
and tail are found both for the proton and neutron, with a positive vector cloud in case of the proton 
and a negative one for the neutron, pointing to first details of the nucleon sub-structure. However, the 
specific nature of these distributions remained yet unexplained. 
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Figure 1: Charge distribution for the proton and the neutron implied by the form factors of Ref. [10 
from [10J, ® (1961) by the American Physical Society. 



Figure 2(b); 



In 1964 Gell-Mann [TT] and Zweig p2] proposed the quarks as building blocks of hadrons to catalog 
the plethora of observed mesons and baryons. During the late 1960ies the MIT-SLAC experiments 
[T3ltT9"] measured deep-inelastic electron-nucleon scattering at the Stanford Linear Accelerator at much 
shorter distances and beyond the resonance region. The important finding of these experiments were 
scaling and the observation that the longitudinal structure function is small, Figure [2j confirming 
a prediction by Callan and Gross [20] for scattering off spin 1/2 particles. The scaling behaviour 
of structure functions had been predicted by Bjorken using current algebra methods [21] ■ These new 
observations led Feynman to the parton model [22U23] of point-like fermionic constituents of the nucleons 
which react at high virtualities with the exchanged gauge bosons in the deep-inelastic process directly. 
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Figure 2: Left: An early observation of scaling: vW% for the proton as a function of — q 2 for W > 2GeV, at 
x = 1/uj = 0.25; Right: The Callan-Gross relation: K$ = F2/{2xF\) — 1 vs — q 2 . These results established the spin 
of the partons as 1/2; from [18], ® (1991) by the American Physical Society. 

Deep-inelastic scattering off constituent quarks has been discussed as early as 1967 [21] in connection to 
data of that time [25] . After the discovery of scaling at SLAC also data taken in other experiments were 
analyzed for this behaviour. One example concerns data taken at DESY at lower values of \q 2 \ [26], cf. 
Figure El presented using the Rittenberg- Rubinstein variable ujw 

1 G. Zweig named the hadron constituents aces. 
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Figure 3: The function i/W 2 plotted vs |g 2 | assuming i? = \q 2 \/v 2 for different fixed uj w = (1/x + M| r /|(/ 2 |)/(1 + 
0.2/Iq 2 !); from [26] ©(1972) by Elsevier Science. 

The parton model introduced a new level of compositeness for fermions being confined inside hadrons 
and related to the strong interactions. The final quantum field theory of the strong interactions devel- 
oped over a series of years. Already in 1965 Nambu [27] proposed a Yang-Mills [28] SU(3) gauge theory 
for the strong interactions, based on a three- valued charge degree of freedom [29]. Before a symmetry 
was introduced using para-statistics [30] which later became color. At this time it was unknown whether 
Yang-Mills theories could be renormalized. The formalism by Faddeev and Popov [31] needed for their 
quantization in covariant gauges has been found two years later only. The renormalization of massless 
Yang-Mills theories was proven by 't Hooft in 1971 [32] and Quantum Chromodynamics (QCD) as 
the theory of strong interactions was proposed by Fritzsch and Gell-Mann in 1972 [33] and Fritzsch, 
Gell-Mann and Leutwyler [31]. In 1973 Gross, Wilczek [35] and Politzer [36J studied the running of the 
strong coupling constant of color octet Yang-Mills theory with color triplet quarks and found asymptotic 
freedom, see also J27JEB]- The Lagrangian of QCD, referring to the covariant i?^-gauges, is given by [2H] 

£qcd = W-™*] - \*r - ^(9^ a (x)f + d, X «(x)D ab >» Xb (x) , (1.2) 

where ip q (x) denotes the quark fields, A%(x) the gluon fields, Fg v = d^A v a — d v A^ a + gf abc A b ^A c > u the 
field strength tensor, f a b c the structure constants of SU(3) C , the gauge group of QCD, £ e R the 
gauge parameter, Xa{%) the ghost field, and the covariant derivatives D ab = <5 a6 c^ — gf abc A c ^(x), p^ k = 
Jti [Sjkd^ 1 — igA^tjf^ , with t a the generators of SU(3) C . Based on this, perturbative calculations in 
Quantum Chromodynamics can be performed at large virtualities. Due to their high complexity these 
calculations are usually being performed using computer algebra programs, a first dedicated of which 
was SCHOONSCHIP by M. Veltman [U]. 

At short distances, resp. large scales of the momentum transfer, the nucleon structure functions 
Fi(x, Q 2 ) obey the light-cone expansion jJTHH] through which they are represented in terms of Wilson 
coefficients C\ and non-perturbative operator matrix elements (OMEs) of local operators, charac- 
terized by their twist [35]. At leading twist the functions fk are the parton densities and at higher 
twist they correspond to parton correlation functions. In this separation an arbitrary scale fi 2 F , the 
factorization scale, emerges. Applying the renormalization group equations [3BT - H9] to the structure 
functions one obtains separate evolution equations for the Wilson coefficients and operator matrix ele- 
ments which are usually of matrix-type due to the mixing of different operators. The solution of these 
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evolution equations exhibits scaling violations □, implied by the renormalization group equations, and 
the running of the strong coupling constant and finite quark masses. 

During the last 40 years the deep-inelastic structure functions both in neutrino and charged lepton- 
nucleon scattering were measured at a steadily increasing precision and allow at present for important 
QCD precision tests at the per cent level. The strong coupling constant a s (M§) can be measured at an 
accuracy of ~ 1% from the scaling violations of the deep-inelastic structure functions. Very accurate 
extractions of the twist- 2 parton densities are possible, which are an essential ingredient for the physics 
at hadron colliders as Tevatron and the LHC, and are thus instrumental in the search for new elementary 
particles. 

At the theoretical side, the deep-inelastic structure functions are best described at the level of the 
twist- 2 contributions. The QCD coupling constant is known to 4-loop order [S21IM]- The anomalous 
dimensions and Wilson coefficients for the unpolarized massless case were calculated to 3-loop order 
[55t - f6T] . There are first results at 4-loop order [62TI64"] . Heavy quark corrections are known to 2- loop 
order [651 EE] and for larger scales a series of Mellin moments has been calculated to 3-loop order 
[67] . In the polarized case the level of 2-loop corrections has been reached for massless and massive 
corrections [68TP72"] . with first results at 3-loop order [73J. The present massless and massive 3- and 4- 
loop calculations require powerful computer algebra systems and packages like FORM |74j, MATAD [75] . and 
Sigma [76], cf. also [77] , to perform these voluminous calculations which amount even in parallel form 
to several CPU years and request 0(100 Gbyte) RAM at multi-processor systems. QED corrections to 
the deep-inelastic scattering cross sections are large in certain kinematic regions and have to be known 
at the accuracy defined by the measured data. The region of small values of x or large hadronic masses 
at a given virtuality requires special attention because of potentially large perturbative corrections. 
In the large x region target mass effects and higher twist corrections contribute. For large values of 
x higher order resummations are important. In the polarized case the structure function g2(x,Q 2 ) 
provides an inclusive observable at which twist-3 contributions can be studied. Structure functions 
and parton distributions obey a series of sum rules, to which also QCD corrections are available and 
one may perform QCD tests using these relations. The nucleon spin problem [78J rouse the question 
of also angular momentum contributions. These can be accessed in deeply-virtual Compton scattering 
and similar reactions [79] as has been shown in [80] . 

In this article we review the present theoretical understanding of the deep-inelastic process, mainly 
in the region of large scales. In Section [2] the scattering cross sections are described and a brief survey on 
the present deep inelastic data is given. Section 3 deals with the light-cone expansion and in Section 4 
the parton model is described. The renormalization of deep-inelastic structure functions in the twist-2 
approximation is discussed in Section 5. Related to this, a survey on the perturbative expansion of the 
strong coupling constant a s (Q 2 ), the anomalous dimensions, resp. splitting functions, and the massless 
Wilson coefficients are given in Section 6-8. In Section 9 we discuss QED and electro- weak corrections to 
deep inelastic scattering. The heavy quark corrections are described in Section 10. Section 11 deals with 
the target mass target mass corrections. The solution of the evolution equations is outlined in Section 12. 
Here we also discuss the results of different NNLO QCD analyses and the status of the determination 
of a s (M|) from deep- inelastic data. In Section 13 a series of aspects of small- a; resummations are 
considered. Resummations in the region of large values of x are discussed in Section [TD Major sum 
rules and integral relations for deep-inelastic structure functions are summarized in Section 15. Higher 
twist corrections and aspects of nuclear parton distributions are considered in Sections 16 and 17. In 
this review we will fully concentrate on firmly established results on the basis of perturbative Quantum 
Chromodynamics and will not discuss other approaches of a more phenomenological character, including 
other assumptions. 

2 Scaling violations emerge naturally through radiative corrections in all renormalizable field theories. Drell and 
collaborators at the end of the 1960ies were seeking for scaling describing structure functions by fermion- meson interactions 
with loop corrections and found scaling violations in general [SDJETJ, see also [52] . 
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2 The deep-inelastic process 



The deep-inelastic process, at Born level, can be illustrated by the diagram shown in Figure 4. A 
lepton I = e ± /i ± , Vi(pi) with momentum k\ scatters off a nucleon N exchanging an electro- weak gauge 
boson V — 7, Z°, to a lepton /' with momentum k 2 = k\ — q and an ensemble of hadrons X in an 
effective 2 — > 2 process k\ + p\ — > k 2 + p 2 . Here pi is the momentum of the incoming nucleon and p 2 
the momentum of the outgoing hadrons. 

In many deep-inelastic scattering experiments the kinematic variables can be measured only from 
a sub-set of the momenta of the external particles. This applies in particular to the momentum of 
outgoing neutrinos but also for the momentum p 2 in fixed target experiments. As will be shown later in 
Section [9] the size of the radiative corrections strongly depends on the choice of the kinematic variables. 
For the measurement at HERA a wide variety of sets of kinematic variables was designed to allow for 
cross checks and to limit the size of QED radiative corrections. 




Figure 4: Diagram describing deep inelastic IN scattering. 

At Born level the unpolarized scattering cross section depends on the virtuality of the exchanged gauge 
boson V, q 2 = —Q 2 and the inelasticity y. At a given cms energy s = {k\ + pi) 2 these quantities define 
the Bjorken variable x, 

P\.(k\ — ko) , 

y = \ 2- 1 

Pi-h 

x = ^ . (2.2) 

sy 

The different sets of variables are : 

i) Leptonic variables q = q\ = k 2 — kx, yi = pi-(ki — k 2 )/pi.ki 

ii) Hadronic variables |M] Q = Qh — P2 — Pi, Vi — Px-{P2 — Pi) / Pi-k\ 
Hi) Jacquet-Blondel variables [82J Q 2 JB = (p 2t j_) 2 /(l — yjB), yjB = S/(2_E(fci)) 

s = J2h( E h-Ph,z) 

iv) Mixed variables [81] q = qi,y m = yjB 

v) Dnvble avale method l8Sl O 2 - AE{k 2 ) 2 cos 2 {6 {k 2 ) / 2) 

v) Double angle method [83] Q DA - + ^{6{k 2 )/2) cos(6(k 2 )/2) tan(% 2 )/2)' 

„ _ 1 sin(g(fca)/2) 

~ sm(6(k 2 )/2) + cos(6(k 2 )/2) tan(^(p 2 )/2) ' 

vi) 9y method [84j Q 2 y = 4E(k 2 ) 2 (l - y JB ) \ ± ^(fcj) - Wv = VJB 

vii) S method [85] Q% = ^L, y s = ^ - E[h)[ f_ co<6[m 

viii) eS method [SS] Q 2 eT , = Q 2 , y e T, = 



Here E(li) and 9(lj) are measured in the detector's rest frame. The invariant mass W of the final state 
hadrons is given by 



Pl = {q + Pi? = M 2 + 2q. Pl - Q 2 = M 2 + Q 2 (l - x)/x 



(2.3) 



and M is the nucleon mass. The inclusive lepton-nucleon process is denoted as deep inelastic if W > 
2 GeV and the values of Q 2 are sufficiently large, normally Q 2 > 4GeV 2 . If W ~ M the process is 
called (quasi-)elastic, allowing also for proton-neutron transitions. Here the scattering cross section is 
described by the (quasi-)elastic nucleon form-factors [86,!H7]. The resonance region is characterized by 
M < W <2QeY . Corresponding models for the scattering cross sections were given in [88J. 
The double differential scattering Born cross section off unpolarized nucleons read 
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F 2 {x, Q 2 ) + Y-xF 3 (x, Q 2 ) + y 2 1 - 



2m 2 l 



2xF 1 (x,Q 2 



I 2 



{Y + W? c (x, Q 2 ) ± yixH/ 3 NC (x, Q 2 ) - y 2 W^ c (x, Q 2 )} 
{Y + W 2 cc (x, Q 2 ) ± Y_xW? c (x, Q 2 ) - y 2 WE c (x, Q 2 )} 



(2.4) 
(2.5) 

(2.6) 



where a and Gp denote the fine-structure and Fermi constants, M^ z and mi are the W, Z-boson and 
lepton masses and 



Y- 



± 



l±{l-y? . 



(2.7) 



Three of the electro- weak parameters in (I2.4H2.6I) are independent. One choice is given by the set 
(a, Mz,Gp). The others are expressed by corresponding relations, cf. [90]. In case of neutral current 
scattering the functions Fj are propagator-weighted structure functions, [90J : 

A, 2 (*,Q 2 ) 
xF 3 (x,Q 2 ) 



F 1>2 (x, Q 2 ) + 2\Q e \(v e + \a e ) X (Q 2 )G lj2 (x, Q 2 ) + 4(v 2 e + a\ + 2\v e a e ) X 2 (Q 2 )H 1 ^x, Q 2 ) 

(2.8) 

-2sign(Q e ) {\Q e \(a e + Xv e ) X (Q 2 )xG 3 (x,Q 2 ) + [2v e a e + X(v 2 + a 2 e )] X 2 (Q 2 )xH 3 (x,Q 2 )} , 

(2.9) 



with Q e = -1,A = sign£,v/ = 1 - A\Q f \ sin 2 9f, a f = l,x{Q 2 ) = G F M 2 Q 2 /{V2 Ma{Q 2 + M 2 )). df 
denotes the effective weak mixing angle, Mz the Z-boson mass, and £ the lepton polarization. The 
structure function 2xF x can be expressed using 



2xF 1 (x, Q 2 ) = F 2 (x, Q 2 ) - F L (x, Q 2 



(2.10) 



and similar for G\,H%,Wi . At this level the Callan-Gross relation [20] implies vanishing longitudi- 
nal structure functions, which take finite values due to QCD corrections and target mass effects, see 
Sections IBlfTTl 

The scattering cross sections fl2.4H2.6p can be represented in the general form 



E 



d 3 



dxdydcj) 



j2p,(s,Q 2 )Lrw, 



I 1 



(2.11) 



where Pi(s,Q 2 ) denotes a propagator term and and Wf 11 ' are the respective contributions to the 
leptonic and hadronic tensor. Here we specified also the azimuthal angle of the final state lepton which 
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emerges in the scattering cross section of polarized leptons off transversely polarized nucleons. The 
leptonic tensor can be calculated perturbatively. The hadronic tensor is given by 



WW = dV*(ps | [j; i (^) t ,^ 2 (o)] I ps) . 



(2.12) 



Here S denotes the four-vector of the nucleon spin with S ■ P — and the normalization S 2 = —M 2 . 
In the framework of the quark-parton model the currents are given by 



/,/' 



(2.13) 



with v J q and a? q the vector and axial- vector couplings of the quarks. For charged current interactions Uff 
denotes the Cabibbo-Kobayashi-Maskawa matrix, whereas for neutral current interactions Uff = Sff. 
The hadronic tensor obeys the following properties [9T|[92] 



Covariance 
Hermiticity 

Spectrality 
Causality 



W^iq', P') = AJ,A£WV(g, P), A g L\ 
WUq,P) = W*(q,P) 



> 1 



W lu ,{q,P) = 0, -, 
W ia ,(x,P)= [ dfiqe-^W^fap) 



T - invariance : W^(q, P, S) = [W^q, P, S)}* 
where L + is the orthochronous Lorentz group and = a? . In general, 



Parity 
Symmetry 
Current conservation 



W tMU (q,P,-S) = W^(q,P,S) 
W tMU (q,P) = W^(-q,P) 
q^W^ = W» v q v = . 



(2.14) 
(2.15) 

(2.16) 

(2.17) 
(2.18) 



(2.19) 
(2.20) 
(2.21) 



are not obeyed. We construct the following hadronic tensor including the case of polarized targets, [93J : 



W, 



9^ + 2 



P P a x P a 



F x (x, Q 2 ) + -g-^Fzix, Q 2 ) - is^xo 



P ■ q 



2P-q 



F 3 (x,Q 2 



^F 4 (x,Q 2 ) + 



(p»Q» + Pvq^) 



p ■ q 



2P-q 



F 5 (x,Q 2 



q x S a 



+i£ fJ ,uXa^ — 9i(x,Q )+ie fluX a 
P ■ q 



+ 



P Q i c p p p 

— — S ■ q- 



2 

P. A 



P ■ q 



q x (P -qS° -S- qP a ] 

(p~W 2 

93(x, Q 2 ) 
P-q 



92(x,Q 2 



(P-q) 2 ' 

PcrS\ 

P ■ q 

v + q^p )S ■ q 



q 2 P-q 

+i£»u\*j^9e{x, Q 2 ) + S ■ 9 ,p Q \ 2 97(x, Q 2 ) 



with 



P = P 



2(P-q) 2 
P-q 



l^WJ + 2P-g 



l 99(x,Q 2 



(2.22) 



Q 1 



-Qt* 



Sfj, — Sfj, 



S-q 



The choice of the polarized structure functions beyond is not unique in the literature, cf. |94j . For 
the functions #3,4,5 we used the convention of [HI], and for the further structure functions that of [HSUHS] ■ 
The spin- vector in the longitudinal and transversal case is defined by 



S L = (0,0, 0,M) , 

St = M(0, cos a, sin a, 0) 



(2.23) 



For the purely polarized nucleon contributions one obtains the following differential scattering cross 
sections, 193 1 : 



d 2 a{X,±S L ) 
dxdy 



d 3 a(X, ±S T ) 
dxdyd<p 



x 



, / 2xyM 2 
-2Xy \2-y- 



2 1 + 



2xM 2 



S 



i-y 



vx 2 M 2 ■ AxM 2 / xyM 2 ' 
xg\ + 8\?—z—g % 2 + — [l-y- JL — j g 3 



xyM 2 



S 



g\ - 2xy z 1 + 



2xM 2 



9~ 



xyM 2 . ( xyM 2 
+ 4A^ gl-2[l-y--» I g"% 




1 / 2xyM 2 \ ■ 2 , 

- - 2 - y - )g\ + - 1 - y - 

y V s J y \ 



cos(a — (ft) [— 2\yxg\ — AXxg l 2 
xyM 2 



g\ + 2xyg\ - 2Xg\ - g\ 



(2.24) 



(2.25) 



Here A denotes the lepton polarization and C 1 



4Qzsin 2 ^ 



■1/2,9$ 



w- 



w - 



1 n n 
L i 9 A 

\z\> 



■hC z = 



= 9v + ^gA,Qv = 1 ,9a = °i9v = 
(g v + Xg A f,C w± = (1±A),^I 2 = 
,,, f z = G F Q 2 /(2V27ra)(M 2 /(Q 2 + M 2 ), V \ Z \ 2 = (rft*)*,^? = [G F Q 2 /{A^a){M 2 w /{Q 2 + M 2 ,)] 2 . | 
The structure functions gj to gg do not contribute in case of vanishing lepton masses and the structure 
functions F 4j5 are related to -^1,2,3- 

The individual structure functions can be measured form the scattering cross sections 02.4H2.6f2.24[ 
12.25ft applying the QED and electro-weak radiative corrections to the data, cf. Section |9j Current 
data analyses in the unpolarized case refer to the structure function data on F2 and charged current 
scattering cross sections for proton and deuteron targets taken at SLAC [98,99], from BCDMS |100pi01j . 
NMC |102l H03j at CERN, and HERA |104j . Surveys on earlier neutrino and charged lepton deep- 
inelastic data were given in |105j . 

Moreover, there are data on the longitudinal structure function Fl(x,Q 2 ) [99 |I100[ [T03, 106-109J. 
To separate the contributions due to the different sea-quark flavors Drell-Yan data [HOj and di-muon 
data [11 lj are used. In some analyses |112[I113] also other sets of inclusive and semi- inclusive deep- 
inelastic scattering data, including data on F 2 QQ (x,Q 2 ), Q = c,b, data from E665, CHORUS, and 
Tevatron data on weak boson production |114] are used. Various PDF-fitting groups use also the 
Tevatron jet data |115] in their analysis. 



3 A factor of 1/2 has to be corrected for 77 



|w±| 2 



i, cf. m\- 
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Figure 5: Left: HERA combined NC e + p reduced cross section as a function of Q 2 for six x-bins compared to 
the separate HI and ZEUS data input to the averaging procedure. The error bars indicate the total experimental 
uncertainty. The individual measurements are displaced horizontally for better visibility; Right : The combined data 
with the HERAPDF1.0 fit |116) is superimposed. The bands represent the total uncertainty of the fit. Dashed lines 
are shown for Q 2 values not included in the QCD analysis; from [104J ®(2009) Springer Verlag. 




Figure 6: The data on Fl versus x obtained by the HI collaboration [109J confronted with the 3-flavor scheme 
NNLO predictions based on the different parton distributions functions (PDFs). Solid line: ABM11 [117| . dashes: 
JR09 [TTB], dots: MSTW [IT2j ). The NLO predictions based on the 3-flavor NN21 PDFs [ffi5] are given for 
comparison (dashed dots). The value of Q 2 for the data points and the curves in the plot rises with x in the range 
of 1.5 -r 45 GeV 2 ; from [TIT] . 



In the polarized case the data stem from SLAC [T2fM23] EMC, SMC, COMPASS [751 [T25HT28] at 
CERN, CLAS [T30HT32] . and HERMES [EEllEg] at HERA, taken at p, d, and n ( 3 He) targets on the 
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polarization asymmetry A\, the ratio g\jF\ or for the structure function g\. Data on the structure 
function g 2 (x,Q 2 ) were taken by the SLAC experiments, SMC and HERMES [12211135] . 
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Figure 7: HERMES results on xg\ and xgf vs x > shown on separate panels, compared to data from SMC [126J, 
E143 H22], E155 [T2"5l[T2"4] and COMPASS [155]. The error bars represent the sum in quadrature of statistical and 
systematic uncertainties. The HERMES data points shown are statistically correlated by unfolding QED radiative 
and detector smearing effects; the statistical uncertainties shown are obtained from only the diagonal elements 
of the covariance matrix. The E143 and E155 data points are correlated due to the method for correcting for 
QED radiation. For the HERMES data the closed (open) symbols represent values derived by selecting events with 
Q 2 > lGeV 2 (Q 2 < lGeV 2 ); from [134] ©(2007) by the American Physical Society. 

The current precision in the measurement of the unpolarized structure function F 2 (x,Q 2 ) is illus- 
trated in Figure [5J showing the combined HERA data together with the world fixed target data. Also 
the individual measurements of HI and ZEUS are compared. The present experimental errors reach 
the 1% level and the scaling violations of F 2 (x, Q 2 ) are clearly visible over a wide range in Q 2 . 

In Figure [6] recent measurements of the structure function Fl(x,Q 2 ) at low values of x by the HI 
experiment are shown together with predictions by different PDF-fits. The structure function Fl(x, Q 2 ) 
can only be measured changing the beam energy in given bins of (x, Q 2 ) to separate it from F 2 (x, Q 2 ) 
via the ^/-dependence of the scattering cross section. A precise measurement requires high luminosities 
and can be expected at the EIC, and at even smaller values of x at a future high energy ep collider. 

Finally, we illustrate the current precision of the measurement of the polarized structure function 
g±(x, Q 2 ) at p and d targets in FigureO If compared to the structure function F 2 (x, Q 2 ) the experimental 
errors are larger since gi(x, Q 2 ) is measured from a polarization asymmetry. The range in Q 2 of these 
measurements is smaller than in case of F 2 (x, Q 2 ) since all measurements were performed in fixed target 
experiments. 

In Figure [8] the kinematic region which has been probed by the different fixed target experiments 
at CERN, SLAC and the HERA experiments is illustrated. For virtualties Q 2 > 4 GeV 2 values of 
x ~ 4 • 1CT 4 are reached. The highest values of Q 2 reached are ~ 20.000 GeV 2 . There is a proposal for 
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an ep experiment in the LHC ring operating at E p = 7 TeV and E e = 140 GeV, which would extend 
the present x and Q 2 ranges by about 1.5 orders of magnitude [136] . 

The kinematic region for polarized data ranges for 0.005 < x < 0.75 and 1 <Q 2 <70 GeV 2 with a 
kinematic correlation between these two variables. At a future facility like the Electron-Ion Collider 
EIC [13T|.I138] this region will be extended significantly. Moreover, the high luminosity available will 
allow precision measurements also for polarized scattering. 




-7 -6 -5 -4 -3 -2 -1 

10 10 10 10 10 10 10 



Figure 8: The kinematic range in x and Q 2 probed by the deep-inelastic scattering experiments at HERA, CERN, 
and SLAC. The region accessible to a possible future e p experiment with E p = 7 TeV and E e = 140 GeV is also 
shown; by courtesy of M. Klein and E. Lobodzinska. 



3 The Light-Cone Expansion 



We consider the hadronic tensor (I2.12p and limit the discussion for brevity to the case of pure photon 
exchange. It is given by the absorptive part of the forward Compton-amplitude T^ v 



WUP,q) = -\mTUP,q) 



7T 



with 



T^(P,q) = i J *ZJ*(P[rjp(Z)J v (0)\P) 
1 / q^qu 



(3.1) 



(3.2) 



-co ^co 



One may integrate f 13 . 2 1) . |44j . which yields 

T, l/ (P,q)=T flI/ (q 2 ,u) = 4vr / d£ / d\£\ \£\ 

For v — > oo the exponentials behave like 



oo jo 



v 2 — q 2 



x(P|T[J M (x)J,(0)]|P) . (3.4) 



i 0p ±i^u 2 -q 2 \€,\ „ p W/(&±l£l) p ±iM\£\x 



(3.5) 
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Due to the Riemann-Lebesgue theorem [139] the dominant contributions to (j3.4p come from 



|£ Q ±|£||<~, KlS^F' (3.6) 

V MX 



with x the Bjorken variable. Thus 



1 



? = t&-?<-^ (3.7) 

and in the Bjorken limit u, Q 2 — > oo the contributions very near to the light cone dominate. 

At very short distances £ 2 ~ the light-cone expansion of the time ordered product of currents has 
the following representation, [4T1 



lim TJ(x)J(0) ~ C l^^ 2 )U ■ ■ •U^'-' MJV (0,/i 2 ) (3.8) 

Here we consider general currents and 0^'"^ N denote local operators, which are finite as £ 2 — > 0, and 
Q^(£ 2 , /z 2 ) are the corresponding Wilson coefficients. Let Do and Dj be the canonical dimensions of 
the operators and currents and N their global spin. The twist r [45J of the operator is given by 



t = D -N . (3.9) 

The Wilson coefficients behave then like 

/ 1 \ -r/2+dj 

££(£V)«(^) • (s.io) 



The local operators of lowest twist are given by 

^,..^(0) = i^S 



(3.H) 



S q;r;til ,..,^) = i lf - 1 8[fajfY lll D„...D lilr il>] (3.12) 
„..(") = 2i N - 2 8Sp[E%JW„...D w _ 1 F%] . (3.13) 

Here the indices q and g refer to quark and gluon field operators, respectively, and A r denotes the 
Gell-Mann matrix of the corresponding light flavor representation; i/j is the quark field, F£ v the gluonic 
field strength tensor in QCD, the covariant derivative, S the symmetry operator for all Lorentz 
indices and Sp the color-trace, where the index a is the color index in the adjoint representation. In 
the unpolarized case jp = 1 and E^ a = 8%, while in the polarized case 7p = 75 and E^ a = \e^ a , 
where e aia2a3ai is the Levi-Civita symbol. The operators (I3.1H43TT3|) still contain trace terms, which 
have to be subtracted. More generally than just by applying the symmetry operation a rigorous twist 
decomposition can be obtained applying the method of |140j . 

Let us consider the case of twist-2 operators in the following. We form the expectation value of 
the operator (13. 8p between nucleon states and perform the Fourier transform of the Wilson coefficients 
through which 

Q 2 \ A fv 2 \ 1 



N 



T{x,Q*) = ^C«^)A N ^j s .3.11. 



is obtained, where 



/ 



^V^-.-U^V) = ^■■■1mC N (^j (3.15) 

(P|O, WJV (0,/i 2 )|P) = P^...P m A N (^-\ . (3.16) 
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The sum in (13.141) runs over even or odd moments starting with a value iVo > depending on the 
crossing relation of the currents, cf. e.g. [93]. Here /i 2 denotes the factorization scale. This scale 
is arbitrary and cancels between the Wilson coefficients and the operator matrix elements since the 
structure functions do not depend on it. Performing a contour integral in x around the singularities of 
( I3.14p one obtains the expressions for the structure functions in Mellin space : 

F 2 , L (N,Q 2 ) = M[F 2 , L (x,Q 2 )](N) = J^C? (^P) A,n (^j ■ (3.17) 

Eq. ( 13 . 1 7j) shows the factorized form of the deep-inelastic structure function. The factorization theorems 
|141] state this form is remaining under higher order corrections. 
Here the Mellin transformation [142] is given by 

M[A(x)](N) = [ dxx N ^A{x) , (3.18) 
Jo 



with the Mellin-convolution 

A(x)®B(x) = I dxi I dx 2 5(x - xix 2 )A(x 1 )B(x 2 ) (3.19) 
Jo Jo 

M[A(x) g> B(x)](N) = M[A{x)}{N)M[B(x)]{N) . (3.20) 

The Wilson coefficients are perturbatively calculable and account for both the massless and massive 
quark contributions. The operator matrix elements A^ are of non-perturbative nature. In case of twist- 
2 they also denote the moments of the parton distributions for massless quarks and gluons, labeled by 
the index i, 

A, N (^)^m^ 2 ) (3.2i) 



P 2 



being discussed in the following Section. 



4 Parton Models 

The SLAC-MIT experiments [T5],[TB] found the strict correlation between the variables v = p.q/Mu and 
Q 2 for large enough values and at the fixed ratio 

" = ^ = x- 
In the Bjorken limit [21] the structure functions are given by 

lim M^(QV) -> Fi(x) (4.2) 

Q 2 ,v— >oo 

lim uW 2 (Q 2 ,u) -> F 2 (x) . (4.3) 

Q 2 ,v— >oo 

In the (naive) parton model by Feynman [22j[23] one assumes that the correlation (14. ip is exact. A 
nucleon at short distances is assumed to constitute of individual partons which are charged and move 
collinear to the nucleon momentum with momentum fractions X{, Yli x i = 1- Let the parton which 
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interacts with the virtual photon carry charge e\ and momentum fraction x\. The hadronic tensor is 
then given by 



W, 



An 
M 



2E 

dx 1 f(x 1 ) — \M\ 2 27r5((x lP + q) 2 -m 2 ) (4.4) 
o ltj \ 



\M\' 



-Ae\ \p^p\ v + PivP'i^ - g^Pi-q] 



(4.5) 



with p[ = pi + q and m the mass of the parton. f(x\) denotes the number density of the struck parton. 
The (^-distribution, which incorporates the correlation (14. Xp in ( 14.4)) . can be written as 

5{{xip + qf - m 2 ) = 5(Q 2 - 2x lP .q) = jr^Sfa - x) . (4.6) 



2Mv 



After performing the integral in ( 14.41) one obtains 

uW 2 (Q 2 ,u) = J^^xfiix), 2MWMQV) = £e?/ 4 (s) 



(4.7) 



Here, we summed over all charged partons and anti-partons contained in the nucleon. In the Bjorken 
limit the structure functions are thus described by 



(4.8) 



The parton distribution function fi(x) agrees with the Mellin inversion of Eq. (I3.2ip . Furthermore, the 
Callan-Gross relation |20] 



F 2 (x) = 2xF 1 (x) 



(4.9) 



holds. The above argument could have been also given assuming a micro-canonical ensemble in ( 14.41) . 
In the parton model the unpolarized structure functions at the Born level in ( l2.4H2.6p are given by : 



F 2 (x,Q 2 ) 
G 2 (x,Q 2 ) 
H 2 (x,Q 2 ) 
xG 3 (x,Q 2 ) 
xH 3 (x,Q 2 ) 
W? c (x,Q 2 ) 
W? c (x,Q 2 ) 



x 



W 2 cc > + (x,Q 2 ) 



W^-(x,Q 2 ) 
*Wf c ' + (x,Q 2 ) 
xW^-(x,Q 2 ) 



YlQ 2 M^,Q 2 ) + q(x,Q")} | 7 | 2 
q 

^|Q 9 Mg(x,Q 2 ) + g(x,Q 2 )] \jZ\ 

i2>" + a2 Mx,Q 2 ) + q(x,Q 2 )] \Z\ 2 

J2\Q q \a q [q(x,Q 2 )-q(x,Q 2 )} \>yZ\ 
i 

EV 9 [g(x,Q 2 )-g(a;,Q 2 )] \Z\ 2 

zZW + ^^q^ + u^Q 2 )] \ z \ 2 

i 

2x S ^\ i a i [q i {x 1 Q 2 ) - q t (x, Q 2 )} \Z\ 2 

i 

2xY^[d t (x,Q 2 )+u l (x,Q 2 )} \W + \ 2 

i 

2xY,[u l {x,Q 2 ) + d i (x,Q 2 )] \W~\ 2 

i 

2xY,[d t (x,Q 2 )-u l (x,Q 2 )] \W + \ 2 

i 

2xY,Mx,Q 2 )-d l {x,Q 2 )] \W~\ 2 



x 



X 



X 



X 



X 



(4.10) 



14 



Examples for the polarized structure functions (12. 24112. 25 p are [5HH44, 145J: 

g l f(x,Q 2 ) = ^Q 2 [Ag(x,Q 2 ) + Ag(x,Q 2 )] 

g[ Zl \x,Q 2 ) = J2l(v 2 + a 2 q )[Aq(x,Q 2 ) + Aq(x,Q 2 )] (4.11) 
q 

g l 5 Zl \x, g 2 ) = YA v « a « t A ^ x ' $ 2 ) ~ A ^ x ' ^ ' etc - 

q 

Here qi(x,Q 2 ) and qi(x,Q 2 ) denote the quark and antiquark distributions of the up (u) and down (d) 
type for four active flavors in the massless limit. Aqi(x,Q 2 ) and Aqi(x,Q 2 ) are the corresponding 
polarized parton densities. Q q denotes the charge of the quark. Furthermore, 

g 4 (x J Q 2 )=2xg 5 (x,Q 2 ) (4.12) 

holds. The remaining structure functions contain also twist-3 contributions [9411145] . cf. Sections [TBI 16. 

The applicability of the parton picture rests on the comparison of two times |146j : T mt - the inter- 
action time of the virtual gauge boson with the hadron and rufe - the life-time of individual partons. 
Let these times be measured in an infinite momentum frame, with p the large longitudinal momentum. 
Applying old fashioned perturbation theory, cf. |147j . they are given by 

T " ~ vron^x) (4 - 13) 

T "™ ~ E,(J3,--E) = E,(*l, + m?)M-M 2 ' (4 ' 14 ' 

These are non-covariant quantities. Note that in this approach energy is not conserved across a vertex, 
while momentum is conserved. Here, go is the energy component of the exchanged gauge boson and 
Ei denote the energies of the individual quantum fluctuations over the hadronic background and E is 
the average energy. The individual momentum fractions are Xi with Yli x i — 1- m i denote the partonic 
masses and M is the nucleon mass. The parton model is applicable if 

Rt = D*H > i . (4.15) 

r int 

This ratio is a covariant quantity. Yet px has to be a large momentum. If ( 14. 15ft is not fulfilled the 
exchanged gauge boson is unable to resolve an individual parton. 

In the massless case m^, M = and assuming that all partons having the same transverse momentum 
one obtains 

R T = ^Zf\ « QH \- X? (4.16) 



2kU 



The last expression in ( I4.16P is obtained assuming two 'essential' parton with X\ — X, X2 — 1 — x. 
Eq. ( I4.16P is valid for Q 2 3> k\ and if x neither becomes large or very small, sufficiently away from the 
elastic and the high energy region. These arguments are qualitative of course. However, close to the 
excluded regions multi-parton and higher twist effects are expected to contribute essentially. 

The momentum fraction z of the nucleon momentum p carried by the single parton depends on the 
mass of the parton before (mi) and and after (mp) the interaction with the intermediate gauge boson 
as well as the nucleon mass M |148j 

_ Q 2 + m 2 F -mj + y/(Q 2 + m%- mj) 2 + AmjQ 2 
2(u + v 7 ^ 2 + M 2 Q 2 ) 
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Important special cases are M — 0, mi = and M — 0, m/ = m^, 



*(M = 0,mj = 0) 
z(M = 0, m./ = mp) 



x 

x 
2 



1 + 



Q 2 



i + i 1 + 



2m^ 



(4.18) 
(4.19) 



Since z < 1 in both cases x is constrained to values smaller than one. For mj = mp = 0, M ^ 0, z 
coincides with the Nachtmann variable £, (11 1.11) . In the fully massless case the momentum fraction is 
given by the Bjorken variable 



z = x 



(4.20) 



In many calculations the so-called collinear parton model is used. Here the partons carry the mo- 
mentum 



p = zp , z G [0, 1] 



(4.21) 



This Ansatz is sufficient in cases where mass effects or transverse degrees of freedom can be safely 
neglected. 

The covariant parton model generalizes (I4.2ip accounting for transverse momentum effects and the 
finite nucleon mass [145[ [149, 150J. In Sudakov variables the parton momentum k in terms of the 
momenta p and q, with p 2 = M 2 , p.k± = q.k± = 0,k±.kj_ = —k\ is given by 



k = xp + 



k 2 + k 2 - x 2 M 2 



2 



Ax 2 M 2 



+ 2x 2 M 2 



5" 



4x 2 M 2 \ 
1 + —qT~ JQ + xp 



+ k± . 



(4.22) 



Here k 2 denotes the off-shellness of the parton. This representation is used to derive the correct integral 
relations among polarized structure functions which are sensitive to transverse degrees of freedom, like 
g 2 {x,Q 2 ), cf. Sections 15,16. 



5 Renormalization and Factorization and Deep-Inelastic 
Structure Functions 

The quark masses and the strong coupling are scale dependent quantities due to renormalization, while 
observables like the structure functions are independent of these arbitrary scales. The scale dependence 
can be expressed by the operator D(fj, 2 ) [47114^] 

IV) := ^ + ajjjjj - ^Wlsjfl (5 ' 1) 



with 



/?(a^ 2 )) = ^-4^ (5-2) 



^<"<" 2 » - -tS^W- (5 - 3) 
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such that 



[D(fi 2 ) + 7Jl + 7j2 - - n AlA ] F^x, Q 2 ) = . (5.4) 

Referring to (13. 17j) one obtains the following two renormalization group equations for the operator 
matrix elements and the Wilson coefficients 

J2[^(^ + lfr- n ^- n ^]M N ^ 2 ) = (5.5) 

3 

^[lD0i a ) < Sy + 7j 1 +7*-'y& NS ]a(^^) = 0. (5.6) 

Here 7^, 7a and 7 j 12 denote the anomalous dimension of external quarks, gluons, and the currents, 
which can be non-zero if the currents are not conserved, cf. e.g. [151j . Here the scale /j 2 = fi 2 F refers to 
the factorization scale. Furthermore, one may introduce different scales for the renormalization of the 
coupling and/or the masses. Often these more general scale choices are used for rough estimates of the 
remaining uncertainties due to higher order corrections. 

In the twist-2 approximation the deep-inelastic structure functions obey the following generic rep- 
resentation : 



( Q 2 

Fi(x, Q 2 ) = ^ — > 

i V 



in,., iin,. .r ) . : /' ( a s ,z, ) . (")•'« ) 



with a s = a s (fij i ),m c = m c (/z|,),m b = m b (yu|,), and fip,R the factorization and renormalization scales 
and the sum runs over the corresponding partonic combinations. Furthermore, one has to specify the 
scheme in which a s , resp. m C) & are defined. Performing the perturbative calculations to a given order 
one may match the /^-dependence and formulate so-called scheme independent evolution equations for 
observables, cf. |152H155j . 

Let us now consider some aspects of the renormalization of the parton distribution functions and 
the Wilson coefficients. In the calculation we will refer to the MS scheme in D = 4 + e dimensions for 
the regularization of all singularities, expanding in the dimensional parameter e. In the calculation of 
operator matrix elements and Wilson coefficients a universal factor 



S £ = exp 



I {j E - ln(47r)) 



occurs in each loop order, which is set to one, unlike in the MS-scheme. Here je denotes the Euler- 
Mascheroni constant. To account for mass effects usually different schemes will be used in intermediary 
steps. The renormalized parton densities /' correspond to strictly massless partons and the coupling 
constant is that of the MS scheme. So-called heavy quark parton densities will be introduced at a 
later stage. All heavy quark effects are contained in the Wilson coefficients Cy. Let us now define the 
massless flavor part of an (inclusive) structure function Fi(x,Q 2 ), F™ assless (x, Q 2 )- It is given by ( 15. 7ft 
in case the Wilson coefficients do not contain any direct or indirect heavy quark mass effects. Note that 
this definition is different from one sometimes used in experiment, where one requests all final state 
fermions being massless. Correspondingly, the massive part of Fi is given by 

^ massive (x, Q 2 ) = Fi{x, Q 2 ) - i^ massless (x, Q 2 ) . (5.9) 

The local twist-2 operators (I3.11143TT3|) . resp. the parton distribution functions, are renormalized by 

Of = ^(^)Of, O? = 4(^)01 (5.10) 
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in the non-singlet and singlet case respectively. The Z-factors account for the ultraviolet singularities 

oo 'in 

2§(a.,e) = ay + 2>i£^ (5.11) 

k=0 n=l 

and the anomalous dimensions given by 

7 ns = z -i,ns (//2)//2 ^_ z ns (//2)) 7 s = Z-^^Zf^) . (5.12) 

The unrenormalized coupling a s = g 2 /(167r 2 ) is related to the renormalized one by 

a s = Z 2 (a s ,e)a s . (5.13) 

The heavy quark masses in the Wilson coefficients are renormalized either in the on-shell scheme or 
the MS scheme, 

(5.14) 

The functions 8mx^{e) have been calculated in |156j . 

At large enough scales Q 2 ^> m 2 the heavy flavor Wilson coefficients factorize into massive operator 
matrix elements and the massless Wilson coefficients [71]. In case of the structure function F 2 (x,Q 2 ) 
this is the case for Q 2 >10 nig, with tuq the heavy quark mass. In this kinematic regime the renor- 
malization of the heavy flavor Wilson coefficients can be traced back to the massless flavor ones and 
that of the scheme-independent massive operator matrix elements. As has been lined out in [ST] charge 
renormalization is somewhat more complicated and it is useful to first perform it in a MOM-scheme 
using the background field method |157] before one transforms to the MS-scheme. The massive OMEs 
contain collinear singularities due to massless sub-graphs, which have to be factorized by the matrices 
Y^Mjj?). Unlike the massless case, where 



m 



,m 



m 



ke/2 



5m k (e) 



r^V) = %G" 2 )> (5-i5) 

this is not the case for the massive OMEs. The renormalization of the massless Wilson coefficients is 
performed by r^- (/^ 2 ) taken in the massless case, cf. [61] to 3-loop order. It is evident that the fragile 
framework of renormalization scarcely allows for any manipulations on 'phenomenological' grounds ad 
hoc. In any case, in this way a new scheme is chosen. 



6 The strong coupling constant 

The strong coupling a s (fi 2 ) = g 2 /(167r 2 ) = a s (fi 2 ) / An is a central parameter in QCD. It is not an 
observable itself but the various hard scattering processes are parameterized by it in perturbation 
theory. Different renormalization schemes may be chosen, cf. [158] . Most commonly a s {n 2 ) is expressed 
in the MS-scheme [159J. All quark flavors are treated as effectively massless. 
The running coupling is obtained as the solution of the equation 

= -P a 2 s - fa\ - fcat - ftaf + . . . (6.1) 
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So far the contributions to the /3-function have been calculated up to 4-loop order in the MS-scheme, 
where the LO [35H38], NLO [JMED], NNLO [IS21IIS3], and N 3 LO [331IM] terms are given by@ : 



0o 
0i 

02 
03 



11 - ■ 

102 - 
2857 



-N 



38 



N 



5033, T 325 Ar9 
18 ; 54 / ' 



/ 149753 



I" 



6 



+ 3564C:- 



1078361 6508 



162 



+ -^Cs ) N f + 



( 50065 6472 , 



162 



Ar2 1093 Ar3 

A^ 2 H A/1 

J 729 J 



(6.2) 
(6.3) 
(6.4) 

(6.5) 



Here we refer to the color coefficients in SU(3) C and Nf denotes the number of active flavors. The 
solution of (16.1 p reads [165J 



a s (Q 2 ) a s (Q 



2 ^ + 



ln + $W(a s (g 2 ) ; A) - $ (n WQo); A)- 



(6.6) 



The superscript n denotes the term at which the expansion of the /3-function in (16. ip was truncated. 
In NNLO one obtains 



0i 



Wo 



In 



;r 



/3 + /Siac + 2 x 2 



+ 



g - 2/3 /3 2 

00^^0200 - 01 



arctan 



0i + 20 2 x 
V40o02-0'i 



Note that 



N f <5: 
N f = 6: 



40002 
40002 



0i >o 

01 < 0. 



(6.7) 



(6i 



This explains the smaller change in a s (/z 2 ) from NLO to NNLO at scales below the top threshold, if 
compared to the logarithmic corrections characteristic for NLO. 

Often (16. ip is solved expanding the approximate polynomial in a s after the separation of variables. 
Let us define L = ln(/i 2 /A 2 ). Then, according to the convention in [159] . one obtains |166] : 

'01 



a s (n 



0oL 0£L 
'01 



^ln(L) + 



0$L* 



1 



0$L 4 



(ln 2 (L) - ln(L) - l) + § 

Po Po. 

Jl02 



g(-ln 3 (L) + -ln 2 (L) + 21n(L)-- 



81 



(6.9) 



Here A = A(Nf) denotes the QCD scale. In applications the scale fi 2 varies over a wider range and may 
pass flavor thresholds. As a convention, one identifies the scale \x at which a new quark flavor becomes 
active by m q , q = c,b,t. This definition implies matching conditions in a s (fi 2 ). At LO and NLO they 
are given by a s (Nf — 1) = a s (Nf), while in higher orders more specific conditions apply [166H168] . For 
4-loop running the 3-loop matching conditions are 

a s (N f - 1) _ 



a s (N f ) 



C 2 a 2 s (N f ) + C 3 (N f )ai(N f ) 



with C 2 = -14/3, C 3 = -340.7289736 + 16.79813197(A^ - 1) for y, = m q . In Section 
summarize the current status of the determination of cn s (M§) from deep-inelastic data. 



(6.10) 
we will 



4 For Quantum Electrodynamics also the coefficient fa = 195067/486 + (800/3)C 3 + (416/3)C 4 - (6800/3)C 5 has been 
calculated in |164) . 



19 



7 Anomalous Dimensions and Splitting Functions 



The anomalous dimensions of local twist-2 operators jij(a s , N) have the representation 

oo 

7^(a s ,iV) = ^a s ( / . 2 ) 7 :f- 1) (iV) , (7.1) 

k=l 

where N denotes the Mellin variable. They are related to the splitting functions in z-space by 

7 f- 1} (iV) = - C dzz^Pf^iz) (7.2) 
Jo 

uniquely |169| . cf. also |170j . z denotes the collinear momentum fraction of the parton compared to its 
emitting particle. 

Here / labels the three flavor non-singlet (NS 1 * 1 , NS V ) and singlet (S) cases, with partonic transitions 
j —7- i. The three flavor non-singlet splitting functions are given by, cf. [59J, 

7? 9 S ' ± ( A = Y qg (N)±lU N ) ( 7 - 3 ) 
ll S ' Y (N) = % q (N) - %,(N) + Nfi^N) - j s qq (N)) (7.4) 

with 



7«*(#) = l^AN) = S^ v qq (N) + ^(N) (7.5) 
7^. (AO = 7^A N ) = S ^U N ) + ^U N ) ■ ( 7 - 6 ) 

1qq(N) contributes for the first time at NLO and the non-zero difference 7q 9 (A^) — 7^(A r ) emerges 
with NNLO. Up to exceptions the anomalous dimensions are different in the unpolarized and polarized 
cases. Note that there are sometimes different conventions being used in the literature parameterizing 
the anomalous dimensions, splitting functions and Wilson coefficients, due to a different normalization 
of the strong coupling constant or using the operator fid/dfi instead of fi 2 d/dfi 2 in the renormalization 
group equation. Fermion-number conservation implies 

75f-(l) = (7.7) 
7£ S ' V (1) = 0. (7.8) 



7.1 Leading Order 



The leading order QCD anomalous dimensions in the unpolarized case have been calculated in |171f fl73j. 
They are given by 



7j?(A0 



7$ (AO 



7$ (AO 



7g>(JV) 



C f 



4S 1 (JV) 



—AT F Nf 
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(7.9) 
(7.10) 
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resp. m z-space 



CM 
Cm 

CM 
CM 



4T f N f [z 2 + (1 - zf] 
4C A 



z 
1 



+ --2 + z-z 2 

z 



+ 2/3 6(l - z) 



Here, S\{N) denotes the simplest of the nested harmonic sums [174J. 



N 



Sb, ai ,...,a n (N) 



k\b\ J a 1 ,...,a n ^), 



(7.13) 
(7.14) 
(7.15) 

(7.16) 
(7.17) 



The structure of the splitting functions Pqq' ,Pqg' ,P<jq 



k=l 

, Pqg\ Pgq^ has been known from QED, [175] . see also |176[ 
1177] . Following [50j[5T] splitting functions for fermion-pseudoscalar and fermion-abelian vector theories 

£=^$, C = ^A^ (7.18) 

were calculated in [178T I181] . For Quantum Chromodynamics they were computed in Refs. [182H184] 
in the space-like resp. in |185] in the time-like case. In |183[I186] also the individual helicity contri- 
butions were given, |187j . The calculation of the QCD splitting functions [183J was of importance to 
extend the naive parton model to the QCD improved parton model. The notion of splitting functions, 
furthermore, forms a more intuitive picture, if compared to the more formal description obtained from 
QCD corrections to the amplitudes in the light cone expansion, and contributed a lot to the detailed 
understanding of the respective processes. 



The following relations apply for z < 1 : 



P q9 (i 
p 99 (i 



Furthermore, the integral relations 



/ dzz[P qq (z) + P gq [ 

Jo 

dzz [2N f P qg (z) + P gg ( 



*) 

*) 
z) 



(7.19) 
(7.20) 
(7.21) 



(7.22) 
(7.23) 



hold. 



In the polarized case the leading order singlet anomalous dimensions were calculated in PSEIQSS] 
using the operator approach and in Ref. [183] for the splitting functions : 
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(7.24) 
(7.25) 

(7.26) 

(7.27) 



3 The foregoing paper |190j was not fully correct. 
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7.2 Next-to-Leading Order 



The splitting functions at NLO were calculated in |191[I192] . For space-like virtualities they read : 
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with 



P2(N) 
P 3 (N) 

Pa(N) 
Pb(N) 



16 + 48iV + 56iV 2 + 52iV 3 + 25iV 4 + 23N 5 + 14iV 6 + 5N 7 + N 8 (7.33) 
144 + 432iV - 152iV 2 - 900N 3 - 275N 4 + 592iV 5 + 834N 6 + 512iV 7 + 109iV 8 (7.34) 
576 + 1488iV + 1088iV 2 - 536iV 3 - 848N 4 + 71 liV 5 + 2075 A^ 6 + 1949N 7 + 937iV 8 



-288AT 9 + 48N 10 



(7.35) 



-144 - 144N + 772iV 2 + 844N 3 - 109N 4 - 392iV 5 + 134N 6 + 268iV 7 + 67iV 8 (7.36) 
-8-8N- ION 2 - 22iV 3 - 3N 4 + QN 5 + 8N 6 + AN 7 + N s . (7.37) 



In the polarized case the anomalous dimensions were obtained in 
A 7 Jf->«(A0 
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(7.41) 



In the above expressions one may identify the leading behaviour at small x given by poles at N — 1 
resp. iV = 0, cf. Section 13. The leading large x behaviour manifests itself by powers of S\{N) oc ln(iV) 
in the limit iV — > oo, cf. Section 14. 



7.3 3-Loop Order 

The unpolarized 3-loop anomalous dimensions have first been calculated for a series of fixed moments 
in [55H58] and then in complete form in [59],[60] . Independent checks for the moments were obtained 
as a by-product of the calculation of the massive operator matrix elements [67] for all color factors 
oc Tp and for the leading Nf dependence [193- 1195] . In Mellin-space the anomalous dimensions may be 
represented by basic harmonic sums up to weight w = 5 |196j . Analytic continuations for these functions 
to complex values of A" were given in [197]. In z-space they are given by harmonic polylogarithms [198J. 
The corresponding expressions are very lengthly and will not be presented here. At a given value of 
a s their size can be compared to the corresponding values at LO and NLO, cf. Figure M As present 
QCD-anlyses show, the seemingly small difference between the curves at NLO and NNLO account for 
a difference of Aa s (M|) w 0.005, which is essential given the experimental accuracy of ~ 1% being 
reached at present, see Section [121 



r-i i i | i i i i | i i i i |-q 




N N 



Figure 9: The perturbative expansion of the diagonal anomalous dimensions j qq (N) and ^gg(N) for four flavors 
at a s = 0.2. The pure-singlet (ps) contribution to j qq (N) is shown separately. From |60J, ©(2004) by Elsevier 
Science. 

First numeric results on the quarkonic contributions to the 3-loop polarized anomalous dimensions were 
given in [73] . 
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7.4 4-Loop Order 



At the 4-loop level at present the 2nd and 3rd moment of the non-singlet anomalous dimension have 
been calculated [52TIM] . Recently the general color structure of 7( 3 )' NS (iV = 2) has been given in 
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(7.42) 



(7.43) 



In numerical studies the 4-loop 
[T99] 



moments were compared to the result of a Pade approximation 



(3), Pade 

7ns 



(N) 



7ns W 



(7.44) 



for A 7 / = 3, as genuine number of massless flavors. For the second moment the Pade approximation 
and the exact result deviate by 20.9 % and for the third moment by 14.9 % only. In Ref. [499] an 
uncertainty of ±100% has been assumed for (17.441) . 



7.5 Anomalous dimensions in the large Nf limit 



The leading and sub- leading coefficients in 1/Nf can be calculated for the anomalous dimensions and 
the /5-function using the method of Refs. [200J to all orders in the coupling constant in the MS-scheme. 
The result of these computations give very useful predictions for the complete diagrammatic calculation, 
which is more difficult. Results were obtained for the fermion mass anomalous dimensions |201H203] , the 
QCD /3-function |204] . the unpolarized anomalous dimensions of composite twist- 2 operators |193[ r205j. 
the flavor non-singlet Wilson coefficient of Fi(x,Q 2 ) |206j . the polarized anomalous dimensions |207j, 
and the anomalous dimensions for transversity [208J. Here, the electron and quark mass anomalous 
dimension [20111202] have been calculated to 0(1/Nj). As one example we show the prediction for the 
4-loop flavor non-singlet anomalous dimension 
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(7.45) 



which has been confirmed by an explicit calculation for N = 2 [62,164] and N = 3 [63j . Another example 
concerns a combination of gluonic anomalous dimensions. At 3-loop order one obtains 



(0) 

~(2) , ~(2)7^ 
! 99 ^ 1 91 ~(0) 
799 



C F T 2 R N}< 



64 



(N 2 + N + 2) 



.Q2 

• °1 



Pi(N) 



3 (N + 2)(N + 1) 2 (N -1)N 2 1 27(N + 2)(N + 1) 4 (N - 1)N 4 



25 



64 10iV 6 + 30iV 5 + 109iV 4 + 168iV 3 + 155iV 2 + 76N + 12 
_ ~9~ (iV + 2)(iV + l) 3 (iV- l)iV3 ' V 
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27 (JV + 2)(N + l) 3 (iV - l)iV 3 
Pl (JV) = 33iV 10 + 165iV 9 - 32iV 8 - 1118iV 7 - 5807iV 6 - 12815iV 5 - 16762iV 4 - 13800 N 3 

-7112iV 2 -2112iV- 288 (7.47) 
p 2 (N) = 87N 8 + 348iV 7 + 848iV 6 + 1326iV 5 + 2609iV 4 + 3414iV 3 + 2632iV 2 + 1088iV + 192 . 

(7.48) 

Here 7^ denotes the respective part of highest power in Nf. The result was confirmed in [60], and more 
recently also in [195J. 

8 Coefficient functions 

The Wilson coefficients in Mellin space obey the following perturbative expansion 

Ci(N,a 8 ) = S iq + ^a s (/i 2 )cf } (N) . (8.1) 

k=l 

Similar expressions are also obtained for other hard scattering inclusive cross sections depending on 
a single momentum-fraction scale like for the Drell-Yan and hadronic Higgs-boson production cross 
section in the heavy top-quark approximation |2091I210| . 

8.1 First Order 

The first order unpolarized coefficient functions in the massless case in the MS-scheme are given by, 
cf. [152] : 

C W (N) ( 9iV 3 + 2iV 2 - 5N - 2 3iV 2 + 3iV-2 , 

r (D (N) N ( -2(-2-iV-4iV 2 + iV 3 ) 2(2 + iV + iV 2 ) c 

Likewise the LO massless polarized Wilson coefficients for the twist-2 structure function gi(x,Q 2 ) 
read, [2TTd214] . see also (HE] : 

ACS(N) = C£>(N) (8.6) 



(8.7) 

The first moment of AC[^(N) vanishes, which is also observed for AC[ 2 g(N), [68] . 
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8.2 Higher Orders 

The massless Wilson coefficients for the unpolarized structure functions at 0(a 2 ) were calculated in 
[19211215] for F 2 (x, Q 2 ), in [r9~2ll2T61l2TT] for F L (x, Q 2 ), and for xF 3 (x, Q 2 ) in [2TS]. A series of moments 
was calculated in [219J. The polarized 0(a 2 ) Wilson coefficients for the structure function gi(x,Q 2 ) 
were computed in [68]. All Wilson coefficients can be expressed in terms of harmonic sums up to weight 
w = 4. As an example we show the coefficient functions in case of F F (x,Q 2 ), with the quarkonic, the 
pure-singlet and the gluonic contribution : 
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2N U + 41N 10 + 226N 9 + 556N 8 + 963N 7 + 2733N 6 + 716(W 5 + 861(W 4 + 1969N 3 
-2748N 2 - 864N - 216 (8.11) 
1075Af 12 + 1439CW 11 + 73464A^ 10 + 16074CW 9 + 35682N 8 - 516984N 7 - 979012N 6 
-627068N 5 + 84099N 4 + 300258N 3 + 119124N 2 - 648N + 7776 (8.12) 

(8.13) 

85N 12 + 113CW 11 + 5472 A^ 10 + 930(W 9 - 13574 A 8 - 85432 N 7 - 149336N 6 
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-114524iV 5 - 11383iV 4 + 44894iV 3 + 22992iV 2 - 4104iV - 432 (8.14) 
P 9 (N) = 26N 9 + 539N 8 + 32UN 7 + 8465iV 6 + 9342iV 5 + 841iV 4 - 5720iV 3 - 2193iV 2 

+2484iV+1404 (8.15) 
P W (N) = N 10 -13N 9 -39N S + 222N 7 + 1132N 6 + 1787N 5 + 913N 4 + 392N 3 + QA5N 2 

-324iV-108 (8.16) 

The other Wilson coefficients have a similar structure, cf. |220[I221] . 

At 0(al) a series of moments for the Wilson coefficients of the structure functions was calculated 
in [55TI55] . The complete expressions were computed in [5T] and combinations of charged current 
structure functions in [222] . In Figure [10] the relative size of the QCD corrections to two Wilson 
coefficients is illustrated comparing the corrections 0(a s ) to 0(a 3 ). In the flavor non-singlet case the 
corrections for the Wilson coefficients at NNLO are larger than in case of the anomalous dimension for 
large values of N, resp. in the large x region. 
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Figure 10: The perturbative expansion of the non-singlet (left) and gluon (right) iV-space coefficient functions for 
F2/X at a s = 0.2, Nf = 4. 0(a s ) dash-dotted lines, O(o 2 ) dashed lines, O(af) full lines. In the non-singlet case 
LO corresponds to a®, From [51], ©(2005) by Elsevier Science. 

The massless Wilson coefficient at 0(a 3 ) can be represented in terms of harmonic sums up to 
weight w = 6, resp. harmonic polylogarithms in z-space. The corresponding expressions are very large. 
Structural simplifications can be obtained applying algebraic [223] and structural relations [196J of these 
quantities. The Wilson coefficients obey difference equations of order ~ 35 and degree ~ 1000, which 
can be found in principle determining somewhat more than 5000 moments, [220] . The corresponding 
difference equations can be solved with the package Sigma [76] . 

9 QED and Electro-weak Radiative Corrections to Deep- 
Inelastic Scattering 

The QED radiative corrections to the deep-inelastic scattering cross sections become rather large in 
part of the kinematic region due to logarithmic terms of 0(a ln(Q 2 /m 2 )), with m\ the charged lepton 
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mass. The corrections in case of e ± N are larger than for fi^N scattering. They have to be known 
precisely since their unfolding is usually performed for the scattering cross section itself prior to the 
QCD analyses of the deep-inelastic structure functions. 

A first dedicated calculation of the QED radiative corrections to deep inelastic e ± N scattering was 
performed by Mo and Tsai |224] and applied in the analysis of the MIT-SLAC experiments, [T5|IT6]. 
Later calculations were carried out in [225,226j for l^N scattering. 

The detailed knowledge of QED and electroweak radiative corrections was of special importance 
also for the measurements of the electroweak parameters in deep inelastic u(v)N scattering. Early 
calculations were carried out in |227j . More recent calculations have been performed in [228] in relation 
to the NuTeV anomaly [229] . The QED bremsstrahlung corrections apply to the lepton lines and 
the incoming quark line, since the corrections to the inclusive hadronic final state vanish according to 
the Kinoshita-Lee-Nauenberg theorem [230] . In [231] the leading logarithmic corrections in 0(a) were 
calculated, suggesting first to absorb the quarkonic QED corrections into the scaling violations of the 
quark distributions. 

With the advent of HERA the radiative corrections were partly recalculated and dedicated calcula- 
tions for deep inelastic neutral and charged current e ± p scattering were carried out by different groups, 
also using a variety techniques, [8 ^ |90 ||226|I232H250] . These approaches include both semi-analytic cal- 
culations [Ml9Tl l2lMl2^1MTll250] and calculations based on Monte Carlo methods [232lf233lf2%8lf2^9] . 
The virtual electro-weak 1-loop corrections to neutral and charged current deep-inelastic scattering 
were calculated in Refs. [2321 1233. 2351 1236] . The inclusive bremsstrahlung corrections are often pre- 
sented integrating over the phase space of the emitted photons or lepton pairs, which can be performed 
analytically [234TM3] . 

Dominant contributions to the QED radiative corrections may be obtained using leading log (LLA) 
techniques [8^ll226l[231ll238fT242|,l244j . This approach, which is based on the factorization of (collinear) 
fermion mass singularities, allows to determine the terms oc a\n(Q 2 /m 2 ) in a straightforward way for 
different settings of the measured kinematic variables. Also higher order terms were calculated within 
this approach [84 11244] . 

The first order terms are described by: 
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where P^\z) denote the leading order QED splitting functions with i,j G {e, 7}. They are obtained 
from the QCD splitting functions given in Sect. [7] setting Ca = 0,Cp = 1,T F = 1. The scale of the 
correction is set by the logarithm 
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The second order corrections Q((aL e ) 2 ) are: 



d 2 a l,2loop 

dxdy 



-L e 
2tt 



* dzP^(z) 1 9(z - z )3(x, y, z) ^ 



d 2 - 



x=x,y=y,S=S 




+ 



a y 
2tt) 



dz 



L 2 A 2 /\z) 



f=i,i f 



3(x,y,z) 



d 2 a° 



dxdy 



x=x,y=y,S=S 



(9.2) 



29 



Initial State Radiation S = zS 


Kinematics 


Q 2 


y 


zo 


3(x,y,z) 


leptonic variables 
mixed variables 
hadronic variables 

JB variables 

double angle method 

61, Vjb 
E method 
eE method 


zQi 
Ql 

Qjb{ 1 -Vjb)z 


z + yi-l 
z 
Vjb 

ih 

z 

Vjb 

z 

Vda 

Vjb 

z 

Vy. 

ZyeT, 


1 — 2/z 

1 - xiyi 

Vjb 
Vh 
Vjb 


Vi 

z + yi-l 
1 

i 

z 

1 - Vjb 


z - Vjb 
z Q DA 
n2 z(z - y JB ) 

^ 1-vjb 
Ql 
zQf 


1 - x JB (l - y JB ) 


Vjb 

XT, 
XT. 


z - yjB 

z 

z - yjB 

1 - Vjb 

i 

1 


Final State Radiation S = S 


leptonic variables 

mixed variables 
E method 

eE method 


Qi 

z 

2 (1-^(1-2)) 
Ve 2 

Ql 

z 


s 

Vjb 
vtz 

l-lte(l-z) 

VeT,Z 2 
(l-y s (l- Z )r 


l-mii-xi) 

Xm 
z 
E,/ 

z 


yi 

z(z + yi - 1) 
i 

z 

1 

z 2 

l + J/ eS (l-2) 
(l-I/ eE (l-«)> 



Table 1: Scaling properties of various sets of kinematic variables for leptonic initial and final state radiation, cf. [90J. 



This notion reproduces the soft photon terms of complete calculations in leptonic variables (cf. 
e.g. [90J). The shifted variables x,y, the threshold z and the Jacobian d depend on the choice of the 
external kinematic variables, see Table [TJ S denotes the cms energy and the hats refer to variables in 
the sub-system. 

Here the different second order splitting kernels are given by 
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denoting double-photon radiation, scattering of a fermion into a fermion by a collinear photon, and 
collinear fermion pair production. Here, is the mass of the produced fermion, Qf its charge, N c (f) = 
3 for quarks, N c (f) = 1 for leptons, respectively, and <g> denotes the Mellin-convolution (I3.19p . The 
soft-photon exponentiation is performed solving the non-singlet evolution equation in the range z — > 1 
analytically, cf. e.g. |177j . Since the terms up to 0(a 2 ) were taken into account in Eq. f !9.2p already the 
corresponding contributions have to be subtracted. One obtains 
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with C = -3 In [1 - (a/3ir)L e ]. 




Figure 11: a) Radiative corrections in leptonic variables in per cent for E e = 26.8 GeV, E p = 820 GeV. Dotted 
lines: 0(a), dashed lines: 0(a 2 ), solid lines: in addition soft photon exponentiation, from [90J, ©(1995) by Elsevier 
Science; b) Radiative corrections in mixed variables in per cent for E e = 26.8 GeV, E p = 820 GeV. The lines left 
from above to below correspond to x = 0.5,0.1,0.001 in consecutive order. Full lines 0(a) + 0(a 2 L 2 ) + 0(a 2 L); 
dashed lines: O(a); Figure by courtesy of H. Kawamura. 



Finally, there is also a fermion conversion term in Q(a 2 L 2 ), P(z,Q 2 ;e — > e + ) = (j-) L 2 Pet (z), 

In Figure [11] we illustrate the size of QED radiative corrections for two types of the measurement 
of the kinematic variables in the case of neutral current deep inelastic scattering. In the first case 
the variables x and Q 2 are measured from the scattered lepton only, while in the second case mixed 
variables, cf. Table [TJ are applied, resulting into corrections of different type. The radiative corrections 
for leptonic variables become very large at small x and high y. Choosing the double-angle method, one 
obtains corrections which behave rather flat in y, cf. [84J. The different methods to calculate the QED 
radiative corrections have been compared in 0(a) for a variety of kinematic measurements and are well 
understood. In 2nd order so far only LLA results are available [84" ] 1244] for the full set of external 
kinematic variables studied by the HERA experiments. Only in the case of mixed variables also the 
0(a 2 L) corrections were calculated |251j . 

In the region in which the virtuality Q 2 of the exchanged photon becomes very small the observed 
final state consists of a p^-balanced electron-photon pair, with little hadronic activity near the beam- 
pipe. This contribution is called Compton peak and has been studied first in [224j . It is given by 

cf. |239] for leptonic variables. The scale setting is the same as being used in [234] referring to the 
kinematic approach in which the hadronic structure is dealt with inclusively. A more refined expression 
than (19. 8 p was derived in Ref. |250] . LLA higher order corrections are easily obtained applying the 
corresponding radiators. The process was described using the parton model in |238j regulating the 
collinear singularity by finite light quark masses. 
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In the initial calculations of the QED corrections to deep-inelastic scattering based on the parton 
model quark lines were treated as fermion lines with mass m q in the on-shell scheme |232y233[l235t [236j . 
This led to logarithmic corrections of 0(ae 2 ln(Q 2 /m 2 )) behaving flat in y and growing with x to 
~ 7 — 12% for x = 0.9. ..0.99. Here e q denotes the light quark's charge. The correct way in treating 
these contributions is to absorb them into the scaling violations of the partons along with the QCD 
corrections. At leading order one obtains the radiator 
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In this way the correction becomes independent of the values of the light quark masses, cf. [231112 39, 243j. 
The overall effect is bounded by < 1.5% at HERA and hard to be resolved [252] . Later numerical studies 
of the resummation of the leading order QED effect were performed in [253J. 

The integral leptonic QED corrections to polarized deep-inelastic l^N scattering, using the approach 
of form factors, have been calculated in |254] , also presenting the leading log results. Here also the effect 
of the exchange of electro-weak gauge bosons was included, which is of relevance at high energy lepton- 
nucleon colliders with polarized targets. Leading log results were given in |255] . Monte Carlo codes for 
the QED corrections in case of polarized deep-inelastic were designed in |256j . 

The well-known resummation of the soft corrections [177j to all orders in 0((aL) k ) has been extended 
including the hard corrections to 0((aL) 3 ) in [257(1258] , to 0((aL) 5 ) in the flavor non-singlet case 
in [25911260 ], in the unpolarized flavor singlet case in [261] , and in the polarized singlet case in 



These corrections are universal and can be applied as well for other scattering cross sections. For 
cms energies up to the TeV range these corrections are sufficient. However, depending on the process 
corrections to next-to-leading log may be needed. 

Another higher order correction is due to small- x resummation in case of QED. These are the 
resummed non-singlet corrections of 0((a log 2 (a;))'), cf. 



10 Heavy Flavor Wilson Coefficients 



The heavy flavor contributions to the deep-inelastic structure functions are rather large in the unpo- 
larized case for smaller values of x. Due to the different scaling violations of the massless and massive 
contributions to the structure functions the exact knowledge of the heavy flavor Wilson coefficients is 
important. At leading order they were calculated in the unpolarized case in [263] and in the polarized 
case in [214i f264. 265]. They are given by : 
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with v = v/I — Am 2 z/ (Q 2 (l — z)) and m the heavy quark mass. The heavy flavor contribution to the 
structure functions at LO is given by 
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and analogously for gf^(x,Q 2 ,m 2 ). In the polarized case the heavy-flavor contributions to gi(x,Q 2 ) 
and g2(x,Q 2 ) show an oscillatory behavior [265] since the first moment of Hg^ gi vanishes and Hg]g 2 
obeys a Wandzura-Wilczek relation |266j . In the unpolarized case the heavy-flavor Wilson coefficients 
have been calculated at NLO in semi-analytic form in [65j|66]|^ A detailed proof of hard-scattering 
factorization including heavy quarks was given in [268J. 

In the region Q 2 3> m 2 the heavy-flavor Wilson coefficients factorize into massive operator matrix 
elements with massless external lines^ and the corresponding massless Wilson coefficients [71] . For 
the structure function F 2 (x,Q 2 ) this representation holds at the 1% level for Q 2 > 10 m 2 [7TJ, while 
for Fl(x,Q 2 ) this is valid only at much higher scales Q 2 >800 m 2 . The results of [71] were confirmed 
in |270j and the 0{a 2 e) corrections were derived in |271j . In the polarized case the NLO heavy flavor 
Wilson coefficients were computed for the asymptotic region only [7T1172]. At LO and NLO also heavy 
flavor corrections for the charged current case have been calculated [272] . 

Because of the high accuracy to which the structure function F 2 (x, Q 2 ) is measured, the heavy flavor 
Wilson coefficients have to be computed to NNLO, like the massless contributions. This is analytically 
possible using the method of massive operator matrix elements in the region Q 2 > 10 m 2 . At 3- loop order 
the five Wilson coefficients LjJ, l£? , L s g i , H^f, H s g i contribute, cf. [67]. In [67] a series of Mellin moments 
has been calculated for the corresponding massive OMEs ranging up to N = 10, 14, depending on 
the process, applying the code MATAD [73]. First 3-loop results for general values of N were derived 
in [194pi95j . These calculations require completely different technologies. In [194J the Wilson coefficients 
L^f and Lp gi were calculated completely. All logarithmic contributions ex \n k (Q 2 /m 2 ), k — 1, 2, 3 in the 
3-loop case have been computed for general values of iV in |273j . The technology to calculate the 
ladder topologies is already available |274j and first results have been obtained for other more involved 
3-loop topologies. Here codes like Sigma and HarmonicSums [75 ]l275j are essential, which were developed 
further to cope with problems of this kind. In the asymptotic region also the 3-loop contributions for the 
heavy flavor Wilson coefficients to Fl(x,Q 2 ) were calculated |276j . Furthermore, a series of moments 
of the massive OMEs contributing to transversity at 0(a^) were computed in [277] . The heavy flavor 
Wilson coefficients at higher order can be expressed for the pole or running mass, which is implemented 
through the renormalization procedure, cf. Sect. EJ up to 3-loop order in [57]; for a phenomenological 
application at NLO see [278] . In the threshold region of heavy-quark pair production one may estimate 
higher order corrections using soft resummations. They have been studied in the unpolarized case 
in [279] and the polarized case in [280J. 

The above calculations were all performed in the fixed-flavor number scheme. As has been shown 
in [281, 282J this description is sufficient throughout the kinematic range at HERA for charm and 
bottom. Still one may consider to introduce a so-called variable flavor number scheme, in particular 
to match the universal contributions on the way a single massive quark is becoming light. Here one 
assumes that the transition always occurs for one heavy quark at the time. One problem in this context 
is that m 2 /ml ~ 1/10 and one cannot consider charm quarks as massless at [i 2 = m 2 .. The matching 
conditions for the flavor non-singlet and singlet parton distributions have been given at NLO in [283J 
and NNLO in [57] and the transition functions in [HT ] [7T |[72ll283[l284] . with first general N results at 
NNLO in [19411195] . Different descriptions to interpolate between the heavy-quark threshold region and 
the asymptotic region are considered in the literature. As an example, in the BMSN-scheme [282, 283J 
the interpolation in F 2 (x, Q 2 ) is given by 

f qq,bmsn^ g 2> Nf + i) = F ? Q (x, Q 2 , N f ) + if°' asymp (x, Q 2 , Nf + 1) — F 2 Q<5 ' asymp (x, Q 2 , N f ) .(10.5) 

Here ]?QQ> as y m P denotes the heavy flavor structure function in the limit of vanishing power corrections. 
Other schemes were proposed in 12851 1287] . The schemes differ w.r.t. the way how fast a massive 



6 A fast numerical implementation in Mcllin-space has been given in |267j . 
7 Massive OMEs with massive external lines have been dealt with in |269) . 
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quark becomes light. One may infer the correct behaviour from precision data. The matching scales 
are usually chosen by \i = m q . However, this scale is usually process dependent and may turn out to 
be much larger as has been pointed out in [288J comparing exact calculations with the interpolative 
description. In any case, the reordering of terms in the interpolation schemes, cf. |282[l283tl285l T287j . 
have to be compatible with the MS scheme. Otherwise an according scheme-transformation has to be 
provided since in the PDF-fits usually MS parton distribution functions are determined and nearly all 
partonic scattering cross sections in physical applications are computed in this scheme. 

Finally we would like to mention that at 3-loop order also diagrams with two different fermion lines 
contribute [289J. Here the moments N = 2,4, 6 were calculated based on the code qexp [290J. These 
contributions are even universal. However, due to the fact the m;, is not very much larger than m c one 
obtains here double-logarithmic contributions, which may not simply be absorbed into either the charm 
or the bottom distribution. 



11 Target Mass Corrections 

At low 4-momentum transfer Q 2 and large values of x the nucleon mass M modifies the scattering 
kinematics and therefore the nucleon structure functions. The corresponding relations have first been 
worked out by Nachtmann [291j in the unpolarized case@ Besides the Bjorken variable x the new 
variable 

2 - / 2 - / , . 

f = — = = — ii.i 

1 + y/l + Ax 2 M 2 /Q 2 l + r 

emerges. Here the mass of the quark in the scattering process is assumed to be the same, TOj = m/. 
For x <C 1 or Q 2 — > oo both variables approach each other. The corrections for the different structure 
functions can be derived using the local operators as in the massless case, cf. Sect. El |293j . In the 
following we discuss the corrections for deep inelastic scattering off unpolarized and polarized targets. 



Let us define the integrals [93,293 
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Here a generic function F(£) has been introduced, which reproduces the corresponding structure function 
in the massless case containing the respective electro- weak couplings and parton content. As an example, 
the function F in case of the unpolarized structure functions Fi and F 2 is the same. The target mass 



3 The same method has been applied to the polarized case in [292 . 
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corrections can now be written in the following compact form for the unpolarized structure functions 
-Pi,2,3 and the polarized structure functions gi v .. 5 5 in the twist-2 approximation : 

(11.4) 
(11.5) 
(11.6) 

* = 1,2 (11.7) 
i = 3,4 (11.8) 

(11.9) 

These relations easily transform into the corresponding integral representations^ For the implemen- 
tation of the target mass corrections in Mellin-space one may refer to contour integral representations 
directly, see [9"3"|I293] . In the unpolarized case they are given by 



77 


,X, ) 




F 2 


%Q 2 ) 




Fz 


%Q 2 ) 


= O^H^) 


9i 


>,Q 2 ) 


= O 9i G 2 (0, 


9i 


>,Q 2 ) 




95 


%Q 2 ) 





N+2j+l 



1 M 2 N 
+ x- 



2 Q 2 (N + 2j)(N + 2j + l) 

(11.10) 

with afc = Jp 1 dzz k ~~ l F(z). In the limit M 2 — y the structure functions F lj3 contain one power in x less 
than F 2 . Note that the definition of the OMEs differ in the literature. The target mass corrections 
for the longitudinal structure function are obtained by [293J 

2xF L = r 2 F 2 - 2xF 1 . (11.13) 

The corresponding relations in the polarized case were given in [93] • For lower values of x < 0.7 and not 
too small values of Q 2 > 1 GeV 2 the infinite series in (jll.10fHl.12p may be approximated using the first 
five terms to reach sufficient precision. In the large-x domain convergence is reached using up to ~ 250 
terms. This requires economic implementations of the OMEs changing during fits to data, for the 
respective complex values of k along the integration contour, which has been given in [295] recently. 

Target mass corrections lead to a violation of the Callan-Gross relation [20]. Furthermore, one may 
derive integral relations in the polarized case, cf. Sect. d5j In the limit x — > 1 the structure functions 
(jll.4f[TL9j) do not vanish unlike the partonic functions F(x). This has led to numerous discussions in 
the literature [296 , 297J . Effects in this region cannot be considered without of a careful study of higher 
twist effects. In [298] besides the usual target mass corrections contributions of the jet-function at large 
values of x were considered to improve the description in this region. However, this approach has to 
be viewed in comparison with others w.r.t. a consequent twist expansion and the fact, that the whole 
hadronic final state Fock-space has to be summed over. 

Let us mention that in [93] also the corresponding relations for the twist-3 contributions to the 
polarized structure functions were derived. The target mass corrections in case of twist-2 and 3, (jll.7p . 
were also derived in [299] . For non- forward and diffractive scattering they have been calculated in 
PSHHEni] and [302]. 



3 We agree with the results of [293] up to the obvious typos in (4.19) and (4.22), cf. also [294 . 
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12 Solution of the Evolution Equations and Parton Distribu- 
tion Functions 

In Mellin-space the singlet evolution equations can be solved analytically. We follow |303] and express 
the evolution equation (15. 5p in terms of a s = a s (Q 2 )/4:7i. The r.h.s. is expanded in the coupling 
constant. The singlet evolution equation reads : 
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with R = P //3o, Rk = [Pk ~ Si=i A-Rfc-«]//3o0 Here (3k, k > denote the expansion coefficients of 
the /3-function and Pk, k > are the singlet matrices of the splitting functions and 
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with E(a s , N) = J2i=i b( a s, N) + q{a s , N)) , G(a s , N) the flavor singlet and gluon distributions. 

One may obtain the evolution equations for the three flavor non-singlet cases by replacing the 
matrices in (112.11) by scalars. In the singlet case the matrices Rk do not commute in general. Firstly, 
the leading order solution is found by 
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with the starting distribution q(a , N) and 
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Here r± denote the leading order eigenvalues. The general solution is given by 
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The solution up to k = 3 was given in [304] . 
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Potential poles in Uk(N) are canceled by those in t/ _1 in (j!2.6p . The perturbative solution in N fc LO is 
obtained expanding ( 112. 6p up to the fcth common power in a s and a keeping L(a s , a , N). 
The different structure functions in iV-space are represented by 

3 

F,(a s , ao, N) = J2 C n!(<>s, N)q™(a s , a , N) + Cf(a s , N)E(a„, a , N) + Cf (a., N)G(a 8 , a , N), (12.8) 

n=l 

where C^' ,s denote the corresponding Wilson coefficients and q^ s are the flavor non-singlet distri- 
butions. In a consistent representation one also expands fl 1 2 . 8 [) in powers of a s ,ao, similar to (112. 6p . 
matching the factorization scales. All these operations can be performed analytically. 

To transform to x-space a single precise contour integral is carried out numerically around the 
singularities of the problem. Usually these are poles on the real axis for Re(z) < c, with a given 
constant c. If small- x resummations are included the singularities may be located in the complex plane, 
cf. Section [13] The contour integral is given by 

1 f°° 

xf(x) = - dzlm[e^x- c f(N = C)}, (12.9) 
7T Jo 

since usually f*(N) = f(N*) is obeyed, |305] , and C = c + ze l<i> . The solution of the evolution equa- 
tions in Mellin space allows for fast and very precise numerical implementations. In other approaches 
the evolution equations are solved in momentum-fraction space. Public codes for the solution of the 
evolution equations are e.g. QCD-Pegasus [306] . Hoppet [307] . QCDNUM [308] . and Openqcdrad [309] . 

The distributions g^ s (ao, N), E(ao, N) and G(ao, N) are non-perturbative quantities, which have to 
be fitted to the world precision data on deep-inelastic scattering and suitable other hard scattering 
processes. Their shape is a prioiri unknown and the corresponding parameterizations, valid in the 
kinematic region to be analyzed, have to be found. One way to determine these shapes is to express 
the distribution functions in terms of orthogonal polynomials, which has been studied in detail in the 
past [31011311] . The Laguerre polynomials L n [ln(l/x)] converge fastest |3llJ. In general the number of 
orthogonal polynomials needed is too large compared to the possible amount of parameters which can 
be fitted from the data. One way out consists in designing shapes in terms of polynomials or related 
functions, which are suggested by the orthogonal polynomial analysis and extend these forms gradually. 
The NNPDF collaboration uses neural network techniques |312j to find the correct shape of the parton 
densities and to estimate their errors. Due to the finite experimental accuracy and the limited amount 
of data points available it seems that the present number of parameters cannot exceed 30-40 in the 
unpolarized case. In the polarized case this number is even lower. The errors of the statistical and 
systematic errors of the parton distribution functions can be determined in x 2 -analyses. 




Figure 12: The MSTW 2008 NLO PDFs at Q 2 = lOGeV 2 and 10 4 GeV 2 ; from [Tl2] . ©(2009) by Springer Verlag. 
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The valence quark distributions u v = u — u and d v = d — d are constrained by the following sum rules 
in case of the proton : 

dxu v (x,Q 2 ) = 2, f dxd v (x, Q 2 ) = 1 . (12.10) 
Jo 

Furthermore, assuming twist-2 dominance, the momentum sum rule reads : 

/ dxx [£(x,Q 2 ) + G(x,Q 2 )] = 1 . (12.11) 
Jo 

In the following we give a brief summary on the status reached for unpolarizecf^l and polarized twist- 
2 parton distribution functions (PDFs). They are obtained by the QCD analysis of the deep inelastic 
structure functions, supplemented by other hard scattering processes. The latter data sets mainly serve 
the purpose to resolve the flavor structure of the quarkonic sea and also might give better constraints 
on the gluon. 




Figure 13: Left : Ratios of the PDFs at NNLO and NLO by CTEQ for Q 2 = 4GeV 2 ; courtesy P. Nadolsky. Right : 
Comparison of the gluon distribution by NNPDF2.1 at LO, NLO and NNLO at Q 2 = 2GeV 2 ; from []T3], ©(2012) 
by Elsevier Science. 



12.1 Unpolarized Nucleons 

Early QCD analyses of deep-inelastic structure functions have been performed starting in the late 
1970ies in LO, and later in NLO, to derive first principal shapes of the parton densities, cf. |315H318j . 
Analyses of the Dortmund group |316[I319] . the MRS group |320j . and the CTEQ (then Morfm-Tung) 
group |321] followed. In the first analyses the charm quark contributions were either treated as massless 
or at LO, since higher order corrections were not yet know. 

At present the differential scattering cross sections of unpolarized deep-inelastic scattering in the 
massless case are known to NNLO and the heavy flavor corrections to NLO. The present highest order 
parton fits are based on this approximation. We note that some of the data contain significant higher 
twist contributions which have to be quantified, or, if possible for the corresponding reaction, cut away. 
Furthermore, target mass corrections have to be applied, since part of the scaling violations in the large 
x and lower Q 2 region are due to these, cf. |199j . 

At present the six collaborations ABM pH], CTEQ [322], HERAPDF [323], JR [HE], MSTW [TT2] . 
and NNPDF |119] perform NNLO analyses, partly using different data sets. In Figure [12] an overview 
on the principal behaviour of the parton distributions as a function of x are given for two typical scales 

11 For extensive recent reviews see |3 13113 14] . 
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of Q 2 by MSTW at NLO. The distributions rise towards small values of x and become steeper with 
growing values of Q 2 . 



HI and ZEUS 




Figure 14: Left : Comparing the la error bands of the dynamical (dyn) and standard (std) NNLO gluon JR 
distributions [118J at small x for various fixed values of Q 2 . Note that Q 2 = 2 GeV 2 is the input scale of the 
standard fit. The central NLO results are taken from |324| with uncertainties comparable to the ones shown for 
NNLO for Q above 2 GeV 2 . For comparison the 'standard' NNLO results of AMP06 [325] are shown as well. 
The results at Q 2 = 2 and 20 GeV 2 have been multiplied by 0.5 and 2, respectively, as indicated in the figure; 
from |118| . ©(2009) by the American Physical Society. Right : The parton distribution functions from HERAPDFl.O, 
xu v ,xd v ,xS = 2x(U + D),xg, at Q 2 = 10 GeV 2 . The gluon and sea distributions are scaled down by a factor 20. 
The experimental, model and parameterization uncertainties are shown separately; from [104], ©(2009) by Elsevier 
Science. 

Furthermore, the heavier flavors, charm and bottom, contribute above the corresponding thresholds. 
At large values of x the up- and down-quark distributions show a different behaviour due to the valence 
quark contributions and the various sea quark contributions are of different size. In Figure [13] recent 
results of the CTEQ and NNPDF analyses are shown. The size of the ratios of the PDFs by CTEQ 
obtained in NNLO and NLO are compared at a low scale for the different partons showing a particular 
sensitivity in the region of small and large values of x. The gluon density is an important distribution 
for many processes at hadron colliders. A comparison is shown for the shapes obtained for G(x, Qq) in 
the NNPDF analysis from LO to NNLO. Here the gluon distribution is becoming lower at higher orders 
in the small x region, with still a wide error band. In Figure [TH the evolution of the light sea-quark 
distributions of the JR NNLO analysis are shown in the small-a; region. Both the dynamical approach 
and the standard fit do agree rather well. The results of the HERAPDF 1.0 NLO analysis, which are 
based on data taken at HERA only, are shown also quantifying the parameterization uncertainty. It 
turns out, that the latter is of the size of the present experimental uncertainty. Figure [T5l shows the la 
error bands obtained in the ABM11 analysis and compares the different massless parton distribution 
to the results found by the latest JR, MSTW, NNPDF analyses at the scale fi = 2 GeV. While the 
valence distributions and the sum of the up and down sea-quarks do agree rather well in all the four 
fits, there are still important differences in the gluon distribution at lower values of x. In particular 
the gluon distribution of MSTW is taking low values in the small x region, and eventually becomes 
negative unlike in all the other analyses. There are also differences in the d s — u s = d — u distribution 
and for the strange quark density at medium values of x. 

In Figure [161 different lattice determinations of the moment f Q dx[x(u v — d v ) — x(u — d)] are compared 
with the corresponding value obtained from PDF fits, see also [326J. Very similar values are obtained 
by the different fitting groups, cf. |117j . There is still a difference between the lattice results and the 
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value from the PDF analyses, despite lower pion masses are used in the present lattice simulations, 
which has to be understood further. A determination of the gluon momentum with lattice methods has 
recently been performed in |327j . For a review on the status of the calculations of PDF-moments with 
lattice methods see [328]. 



u.=2 GeV, n,=4 




Figure 15: The la error band for the 4-flavor NNLO ABM11 PDFs at the scale of \i = 2 GeV versus x obtained 
in [117| (shaded area) compared with the ones obtained by other groups. Solid lines: JR09 [118J, dashed dots: 
MSTW [112], dashes: NN21 [U9]; from [TIT] . 



A H QCDSFN, = 2 

X ° a RBC/UKQCD N f = 2 + 1 

V ♦ LHPCN, = 2 + 1 

ETMCN, = 2 

• ETMCN, = 2+1+1 

• update ETMCN, = 2 + 1 + 1 

• ABM11 

' old aSj oTfo oXs oSo oSs aio 
ml [GeV 2 ] 



Figure 16: Comparison of lattice computations for the second moment of the non-singlet distribution as a function 
of the pion mass with the result of ABM11 |117| along with the uncertainties of the respective measurement; 
from [117J, by courtesy of V. Drach. 
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Figure 17: Left : The NLO polarized parton BB distributions [529] at the input scale Ql = 4.0 GeV 2 (solid line) 
compared to results obtained by GRSV [3"3"0] (dashed-dotted line), DSSV (long dashed-dotted line) [551] , AAC 
(dashed line) [332J , and LSS (long dashed line) [333J. The shaded areas represent the fully correlated la error 
bands calculated by Gaussian error propagation; from [329J, ®(2010) by Elsevier Science. Right : The polarized 
DSSV PDFs [331] of the proton at Q 2 = 10 GeV 2 in the MS scheme, along with their Ax 2 = 1 uncertainty bands 
computed with Lagrange multipliers and the improved Hessian approach; from [331 j . ®(2009) by the American 
Physical Society. 

12.2 Polarized Nucleons 

Early determinations of the polarized parton densities were carried out at the beginning of the 1990ies 
at LO under a series of model assumptions, see e.g. [3341 1335] . With the advent of the NLO polarized 
anomalous dimension NLO analyses were carried out |215p333p334p336[l337j ; cf. also |338(.I339] . In these 
analyses SU(3)f symmetry was assumed for the sea quarks. In FigurefTTIthe results of different polarized 
PDF-analyses are compared (left panel). The errors on these distributions are larger than in the 
unpolarized case, with a basic agreement in the results for the valence distributions. For the sea-quark 
and gluon distribution the errors are larger and there is a stronger variation of the current predictions, 
widely within the present errors. We note also other approaches based on statistical distributions [340J, 
which allow for a very efficient modeling. 

To resolve the flavor dependence of the sea-quarks data on semi- inclusive measurements, as e.g. [341J, 
are used. Furthermore, data on photo- and electro-production of hadrons and charm and proton-proton 
collisions at RHIC can be used, along with the structure function data. A recent analysis based on 
these data was carried out in |331j . The semi-inclusive data were also used in the analysis |342j . The 
polarized PDFs of [331J are shown in Figure [T7] (right panel) and allow, yet with larger errors, to derive 
different shapes for the polarized sea quarks. The gluon distribution is found to be much lower than in 
the case of the analysis of only the deep-inelastic structure function data. 

12.3 a a (M|) 

The strong coupling constant can be determined as one parameter in the QCD-analysis of the deep 
inelastic world data along with the parameters of the non-perturbative parton distributions at a given 
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scale Qq. In the following we will compare the results on cn s (M§) in different NNLO analyses. The 
results are summarized in Tabled cf. |117j . 





a s {ivi z ) 
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n i i oa + 0.0019 
_ 0021 
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valence analysis, NNLO |3 1 3 
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U.lloO ± U.UU14 


WO- T?T?1\TQ r, Q lOQOl 
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ArjivlVl 


n 1 1 on _i_ n nm A 
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n(ai. rjoiviiN-approacn [zozj 
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U.11Z4 ± U.UUZU 


dvnamical approach |118J 
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cfanrlaT-rl fif 11181 

staiiuaiu nt [J-J-o] 


ABM11 


0.1134 ±0.0011 


11171 

11 


MSTW 


0.1171 ±0.0014 


13441 


NN21 


0.1173 ±0.0007 


m 


CT10 


0.118 ±0.005 


[346] 


BBG 


n 1 141 + °- 0020 

U.1141 - 0.0022 


valence analysis, 1N°LU(*) [199] 


world average 


0.1183 ±0.0010 


[347] (2011) 



Table 2: Summary of recent NNLO QCD analyses of the DIS world data. 



There are two valence analyses [199(1343] . effectively limited to the region of x>0.3. In [199J an effective 
N 3 LO fit has been performed, including the 3-loop Wilson coefficient, noting that the effect of the 4-loop 
non-singlet anomalous dimension is very small. The other analyses cover the whole kinematic range. 
In the valence analyses and [1171111811282] very similar values of a s (M|) are obtained, which agree 
within their errors. The way in which the heavy flavor corrections are treated implies a systematic error 
of 0.0006. The N 3 LO value is well compatible with the results obtained at NNLO. These values are 
lower than the current world average [347] . Larger values are found in the MSTW [344] and NNPDF 
analyses [345] . A preliminary value has been reported by CTEQ [346J at NNLO, yet with a rather 
large error. In Ref. [117] a detailed analysis has been performed comparing the above results also with 
respect to the pulls given by the different data sets used. Despite the final a s (M|) values of MSTW 
and NNPDF are quite similar, there are still significant differences in the pulls. The Tevatron jet data 
do not cause the larger values. It rather seems that the response of the MSTW and NNPDF fits to the 
SLAC and also the HERA data is causing this difference, which needs to be investigated further. 

In different other reactions larger values of a s (M§) were found. This applies to the analysis of the 
3-jet rate with a s {M\) = 0.1175±0.0025 and inclusive Z-decay yielding a s {M\) = 0.1189±0.0026 
[349,350j. This also applies to the a s (M|) values measured for hadronic r-decays, cf. [347] . On the other 
hand, the analysis of thrust in e + e~ annihilation led to a s (M|) = 0.1153 ± 0.0017 ± 0.0023 [351] , resp., 
a s {M\) = 0.1135 ± 0.0011 ± 0.0006 [3S2]. In a NLO analysis of e+e" ->• 5 jets a s {M\) = 0.1156 ± °Zl 
was obtained. There are also first results on the measurements of a s (M|) with including dynamical 
fermions from the lattice [353]. Very recently a NLO analysis of the ATLAS jet data resulted in 
a s (M|) ~ 0.1156 with a larger error [354J. Similar results are obtained in a multi-jet analysis of e + e~ 
data |355j . Usually the values of a s (M|) obtained at NNLO are lower than those at NLO. In various 
classes of high energy reactions partly larger and lower values of a s (M|) are obtained at present. This 
is due to systematics sources, which have to be understood better in the future. Further clarification 
can be obtained from better hard scattering data. Here the LHC jet data will be of importance in the 
near future. 
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13 Small x Resummations 



Approaching very small values of a; at a fixed virtuality Q 2 the criteria for the parton model, Section HI 
are no longer valid and new phenomena are expected to contribute. The small x region was probed at 
HERA to values of x ~ 10" 4 at Q 2 = 10 GeV 2 and at LHC values of x ~ 1CT 6 at Q 2 = 10 2 GeV 2 can 
be reached. The potentially new effects will influence the size and scaling violations of the structure 
functions both in the unpolarized and polarized case. To which extent these phenomena can be dealt 
with using perturbative methods and resummations based on them, is not finally clear at present. Two 
main extensions of the perturbative approach based on fixed order perturbative QCD were proposed: 

i) the linear BFKL small- x resummation |356] of large logarithms (a s (N c /ir)) k (l/x) \n k (l/x)/k\, and 

ii) saturation models with highly non-linear gluo-dynamics, cf. [357] related to Glauber- models [358] ^ 
In both approaches elements known from the parton picture are used. Since perturbative kernels are 
evaluated it has to be clarified whether factorization holds. In case of the saturation corrections this 
also applies to multi-parton states. Due to the connection between the non-perturbative distributions 
and the perturbative kernels, usually in terms of a convolution, both effects at small and larger values 
of x will contribute. In the following we will not cover the saturation models and related theoretical 
developments in nucleon-nucleus and nucleus-nucleus scattering, being an interesting broad topic in its 
own right. We refer to recent detailed surveys, as e.g. given in [13811361] . and will discuss the effects 
due to BFKL-type resummations. 

The resummation of the leading small x terms [356] are often confronted to those in ln(Q 2 ), stating 
that the first ones results from the strong ordering x% 3> 3> Xi + \... and the second stem from 
strong ordering in the transverse momentum k± ± <C ■ ■■k± i i <C k±^ + i... along a ladder^ Note that this 
is a gauge-dependent statement and requires to use effective vertices in case of the first approach. The 
second case refers to fixed-order perturbation theory using the renormalization group, as outlined in the 
previous sections. In the first approach one would calculate in a more systematic way the scale-invariant 
contributions to the evolution equation (15. 5p in the massless case setting (3 = 0, i.e. considering the 
strong coupling as const. One obtains 

E^ 2 )=E n M)(^) . (13.1) 

The scale invariant part of the anomalous dimension has the representation 

oo 

„(')„! 



7o fc -n T * = 5>oX (13-2) 



and exponentiates to all orders. The representation (113. ip applies also for higher order resummations 
under the above requirements. To treat running coupling effects one has to systematically account 
for scale-invariance breaking effects. Furthermore, to form observables, the small- x resummed Wilson 
coefficients have to be calculated as well. 

Both in the flavor non-singlet and singlet cases for scattering off unpolarized and polarized targets 
the leading order small-x resummations were derived in form of evolution kernels. They seem to 
correspond to a resummation of the leading order anomalous dimensions, as the comparison in the 
known orders show, because the corresponding singularity in the Wilson coefficients turns out to be one 
order lower |59 |, [60 | [363j . 



12 For early estimates of the transition line between the perturbative and non-perturbative domain and the onset the 
possible of saturation effects see |359|I360) . 

13 One may study this process under a more general point of view and consider angularly ordered emissions covering 
both the above cases [362], which allows for interesting applications through Monte Carlo studies. Whereas this unified 
treatment is possible at LO, higher order corrections cannot be cast into this form in general. 
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13.1 The Non-singlet and Polarized Singlet Contribution 



The most singular contributions to the Mellin transforms of the structure-function evolution kernels at 
all orders in a stem from the poles at iV = in each individual order. In the flavor non-singlet case 
the kernels can be obtained from the positive and negative signature amplitudes studied in [364J 
for QCD via 

M [KtM] (iV) = j dxx^K^^a) = ~Tt*{N,a), (13.3) 

with 

^ (N,a) = -nI 1 -]^^-} (13-4) 



x->0 



(N,a) 



-N < 1 



8aC F 



(13.5) 



where z = N/ ^2N c a, D p (z) denotes the function of the parabolic cylinder, and N c = 3 in case of QCD. 

The LO small x evolution kernels in the case of the polarized singlet evolution were derived in [365J. 
The resummed splitting function is given by 



P(x,a s ) = yp 



(0 



1=0 



x-fO a s 



l+l 



log 21 X 



rM-^FoC/V.tOKo;). 



(13.6) 



The matrix valued function Fq(N, a s ) is obtained as the solution of 



F {N, a s 
with 

F 8 (N, a s 



C F 
2C F 



16tt 2 ^M 8 



-2T R N f 
4C A 



Cf — Ca/2 
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-T R N f 
2C A 



2a a 
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d 



-^—-F 2 Q (N, a s 



8ir 2 N 



;i3.7) 



1 1 



+ W C *dN F ^ a ° ) + 8^N 



F 2 s (N,a s ). 
(13.8) 



Eq. ffTBTB"]) obeys [555] 

Pg/{T R N f ) = -Pgj/C F (13.9) 

to all orders, where T R = 1/2 and Nf denotes the number of flavors. The leading contributions of 

the fixed order results in LO and NLO (MS) are correctly described. In the supersymmetric limit 
C A = C F = N f = 1,T R = 1/2 the relations 

PS+ P S = P S+ P S (13-10) 



91 



19 



99 



are obeyed for all I and Eq. d 1 3 . 6 [) can be given in a simple analytic form |366j . It is evident from 
Eqs. f ll3.4ffi~3~l)]) that the poles at iV = being present in the individual orders are resummed into 
branch cuts, which usually exhibit a milder singularity. 

To perform numerical studies it is interesting to consider the effect of less singular terms to estimate 
the stability of the leading order resummed terms. In case of the fixed-order anomalous dimensions 
they are known and one may try Ansatze like 



(a) F{N, a s ) ->T{N, a.) - T(l, a s ) (b) T{N, a s )^T{N, a s ){l - N) 

(c) T(N,a s )^ T(N,a s )(l-N) 2 (d) T(N, a s )^T(N, a s )(l - 2N + N 3 ) 



13.11) 
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cf. p03j. It turns out that the corrections to the flavor non-singlet contributions are numerically very 
small |363] . Larger corrections are obtained in the polarized singlet case, see Figure [TSJ However, here 
sub-leading terms e.g. of the type (113.111 (&)) do nearly completely cancel these contributions again. 
Comparable results were obtained in [367] . Large effects as anticipated in |368] are not confirmed. 




Figure 18: The Q 2 evolution of the polarized quark singlet and gluon momentum distributions evolving from 
Ql = 4GeV 2 ; from [366], ©(1996) by Elsevier Science. 
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Figure 19: Real and imaginary part of the LO BFKL anomalous dimension; from |303j . ©(1998) by the American 
Physical Society. 



13.2 Unpolarized Singlet Distributions 

The LO resummation for the evolution kernel of the unpolarized singlet distributions was derived 
in |356j . It was shown in [369] that the eigenvalue 



(N-l) 



-Xohi 



orJV„ 



[2^(1) - HJl) ~ ^(1 - 



'13.12) 



7T 7T 

represents the LO resummed gluon-gluon anomalous dimension 7^ = jgg(N,a 8 ). The resummed LO 
gluon-quark anomalous dimension is given by ^^(N, a s ) = (Cf /C 'a)1l and the quarkonic terms do not 
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contribute in 0[(a s / (N — I)) 1 ]- Eq. (113. 12[) can be solved iteratively demanding jl(N, a s ) — > a s j (N — 1) 
as \N\ — > oo for iV G C, which selects the physical branch of the resummed anomalous dimension, 

1l = J 99 ,o(N,a s ) = jf^j |l + 2^( 2l+ll f g + 1 (N,a s )^ . (13.13) 

Here we rewrite a s = N c a s /n. 7l has the serial representation 

— / — \ 4 / — \ 6 / — \ 7 

a s . n/ . ( ot s \ / a s \ 2 / a s 



^>*') = W±l + ^\w±l) +^{n±t) +1 ^{n±t) +"' < 13 ' 14 ' 

Under the above conditions one may calculate 7z(iV, a s ) in the whole complex plane. It is a bounded 
function of p = (N — l)/a s , the singularities of which are branch points |370[I371] at 

P! = 41og2, p 2i3 = -1.41048 ± 1.97212 i. (13.15) 

Its analytic structure is shown in Figure [19j Note, that the resummed form of 7i(iV, a s ) removes all 
the fixed-order pole singularities of Eq. (113.141) into branch cuts. Since the known NLO resummed 
anomalous dimensions are functions of 7i(iV, a s ) which introduce no further singularities the contour 
integral around the singularities of the problem has to cover the three BFKL branch points, the singu- 
larities of the input distributions along the real axis left of N = 1, and the remaining singularities of 
the fixed order anomalous dimensions at the non-positive integers [303,371j. Any finite correction to 
jl may thus lead to essential changes of the corresponding numerical results. Early numerical studies 
on the impact of the LO resummed anomalous dimensions were performed in [372J and more recently 
in Refs. [303ll370ll373] . 

The next-to-leading order resummed anomalous dimensions are given by 




Iqg ~ T. a sJ-F 



1 </.</ 

NL 
99 
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with Tp = TrNj. The quarkonic contributions were calculated in Ref. [374 ], as well as the resummed 
coefficient functions C2(N, a s ) and CL(N,a s ). In |375j cl(N, a s ) has also been calculated, giving a result 
which differs form the one in |374] starting with 0(af). 7^ L was computed in |376[I377] . In the 
DIS-scheme ■y^ fL is found to be an analytic, scale-independent function of 7z,(iV, a s ) and reads 



« s 2 + 3 7 l-3 7 ! [B(l- lh l + lL )f 

yi\,tx s ) = if 

where -R(7) is given by 



NL,DiSr N ) = r n - 1 <»- 1 R( lL ), (13.17) 

qg 6tt 3-2 7l 5(2 + 2 7l ,2-2 7l ) ' 



R W = I ~i~T~\ GXP 9 / dz v M + X °( Z) ■ (13 - 18) 

in [374] . In a re-analysis in |375] 

= 2 n S exp - / dz VK) * VK) + Xo(z) . 13.19 

has been obtained, with a different pre-factor. A future calculation of the 4th order Wilson coefficient 
for the longitudinal structure function in the MS scheme may clarify this question further. 
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The NLO resummed gluon anomalous dimension 7^ L was calculated in the Qo~scheme |378j . Here, 
a scale Ql 3> Aq CD is introduced suppressing k± effects for k 2 < Qq, see also |379] J^I One has to solve 
the Bethe-Salpeter equation 

(N - l)G N (q x , q 2 ) = 5 D - 2 ( qi -q 2 ) + J d D ~ 2 q 3 K( qil q 2 )G N (q 3 , q 2 ) (13.20) 

with 

K{qi, q 2 ) = b D ~ 2 {qi ~ ?z)2w(gi) + K Ieal (qi, q 2 ) + ^virtual (13.21) 
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Figure 20: Different contributions to the resummed splitting function xP gg {x,a s ) in the DIS-scheme; Lx: LO 
BFKL; NLx 9f : NLO BFKL quarkonic contr.; NLx: NLO BFKL; from [350], ©(1998) by World Scientific. 



For q\ 3> q\ one diagonalizes as in the LO case using formally the same Ansatz 



d D - 2 dq 2 K( qi ,q 2 ) (q 2 2 ) 
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Here the scale-invariant LO eigenvalue ct s Xo{l) is supplemented by the NLO correction term 
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13.22) 



(13.23) 



(13.24) 



'For a transformation into the DIS-scheme cf. |303j . 



47 



I — — — — — — 

10° 10 1 10 2 10 3 10 4 10 5 10 6 

1/x 



Figure 21: The gluon splitting functions xP gg plotted with a s = 0.2 and Nf = 4 The curves are (from top to 
bottom for xP gg at x ~ 0.2): fixed order perturbation theory LO (black dashed), NLO (black solid), NNLO (green), 
resummed LO (red dashed) and NLO in QrjMS scheme (red solid) and in the MS scheme (blue); from [381], ©(2008) 
by Elsevier Science. 

Whereas the terms in the first two lines of Eq. (113.231) do contain contributions to the anomalous 
dimension up to 0{a 2 s ) the third line contributes only with three-loop order. The agreement with the 
perturbative results to 2- loop order has been known in the 1990ies. The 0(a? s ) from (I13.23[) . given in 
the DIS-scheme, agree with the perturbative result [60] after transforming it to the MS scheme. The 
test of even further agreement in higher orders is left for the future. Note that the addend —6(^3, being 
numerically large, contains contributions of the gluonic contribution to the trajectory function uj(qf). 
The result given in Ref. |382] was confirmed in a different calculation by Ref. [383J. A departing value 
was reported in |384] 

Numerical results on the impact of the leading and next-to-leading anomalous dimensions and 
coefficient functions were provided in a series of detailed studies, see e.g. |303P370L 380 .385] and references 
therein. The matrix formalism for the solution of the all order evolution equations, extending a first 
approach in Ref. |304] to all orders, both for hadronic and photon structure functions, is described in 
Ref. |303] in detail. The quarkonic contributions lead to a strong enhancement of both F2(x,Q 2 ) and 
Fl(x,Q 2 ) at small x during the evolution. However, already simple choices for the yet unknown less 
singular contributions, cf. (113. lip , diminish these effects sizably so that a final conclusion cannot be 
drawn at present. In the case of the resummed gluon anomalous dimension the NLO contributions 
are found to be extremely large and negative. The large rise due to the LO BFKL term is already 
canceled to the level of the fixed order contributions by the purely quarkonic contribution to 7^ L , see 
Figure [20j Adding also the gluonic contribution leads to negative values for the resummed splitting 
function already for a s = 0.2 and x ~ 0.01 which has to be regarded as unphysical. The LO and 
NLO resummed contributions to the gluon anomalous dimension seem to represent the first terms of a 
diverging series, which might be eventually resummed. 

An interesting approach in this direction was performed in Refs. [381 , 386J , and related work in 
[387, 388J. In [381,386j the coupled evolution equations 

j t f + (N,t) = 1 + (a s (t),N)f + (N,t), ^f + (x,M)= X (d s ,M)f + (x,M), (13.25) 

with t = log(Q 2 /Q§),£ = ln(l/x), f + (x,M) = dtexp(-Mt)f+(x,t), and N the Mellin variable. 
7 + (a s ,AQ is the eigenvalue of the singlet-evolution kernel containing poles at N — 1 and x{^s,M) 



15 Despite the fact that in Ref. |384j the quark contribution to (*j(qf) agrees with [38211383] . it may still be that the 
calculation in 384 addresses a different quantity. 
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denotes the BFKL-kernel, where a s is the operator a s (t — > —d/dM). Eqs. f ll3.25|) give rise to the 
so-called duality relation 

X (a s , 7 + (a s , N)) = N, 7 + (a s , xK M)) = M (13.26) 

in the small- x region. These relations are used to organize the resummation of the small-x contributions, 
which leads to a gradual stabilization of the dominant splitting functions at low x. For the splitting 
function xP gg (x) this is illustrated in Figure [211 including all the presently known information. In the 
MS scheme significantly larger values are found than in the 3-loop fixed order calculation at x ~ 10~ 4 
and below. 

Finally, we would like to mention that the effect of potential sub-leading contributions to the LO 
resummed anomalous dimension may be studied within $ 3 theory in D = 6 dimensions, as a simple 
model with 3-boson interactions at the perturbative level. Here the leading order resummed anomalous 
dimension can be calculated for all values of x solving a Bethe-Salpeter equation |389] . The correspond- 
ing singularity at fixed orders is located at N = — 1 due to the scalar field. In |39Uj the LO small- x 
resummed terms were compared to the complete ladder solution adding the NLO corrections. In this 
case it turns out that the pure small- £ resummed terms do not give the correct result. To further 
consolidate the knowledge on parton distributions in the small-x region it seems to be necessary to 
calculate at least one further resummed series, i.e. NNLO BFKL. 



14 Resummation at Large Values of x 

The dominant contributions to the splitting functions and Wilson coefficients in the large- x region result 
form the terms 



ln fe (l-x) 
1 — x 



5(1 -x), ln'(l-x), k,l e N,k > 0,1 > . 



(14.1) 



This is best seen in Mellin space, where 



k\' 



-M 



ln fe (l - x) 
1 — x 



(N) 



:-i) fc+1 Si,...,i(iv-i; 

k + 1 



M [5(1 — x)) (N) = 1 



^M[ln fc (l-*)](iV) = (-l) fc ^i,...,i(iV) . 



k\ 



(14.2) 

(14.3) 
(14.4) 



The harmonic sums with equal index form a polynomial in single harmonic sums, cf. |174[ [223j. 



k 

Si,...,i(N) = \Y j S 1 ,..., 1 (N)S 1 (N) . 



;i4.s) 



1=1 



k - I 



In the limit \N\ — > oo one obtains 



Si,...,i(N) oc ^f(JV) + 0(S k 1 - 2 (N)S 2 (N)) 



;i4.6) 



and 



Sx(N) ocln(iV) + ^-+ 1 



1 



2N 12N 2 120A^ 4 VA^ 6 



(14.7) 
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with N = N + 'j E . 

Let us consider the splitting functions P, t k and Wilson coefficients C Pji 



Pik(x,a s ) 



J2 a 

3=0 



j+ipU) 



ik 



C Pti (x, a s ) = 5 Pi2(3) 5(l - x) + ^2 a l c p}( x ) ■ 



14.8) 



;i4.9) 



The large x structure for the spitting functions Pik and Wilson coefficients known up to O(a^) [59ti6Tj 



is 



CaAX 



A+l(l 
J'=0 

(2C f )' 
(/-!)! 



x)+ + - x) + C, +1 ln(l - ar) + 0[(1 - x)*- 1 ln m (l - x)} (14.10) 

x)+0(l) (14.11) 



Mln*-'(1 



p ?5 (a;)ln 2 ^ 1 (l-x) + 0[(l-x 



)fc >-l ln 2^2 



'(1-*)], 



14.12) 



cf. |391[ [392j. Here c a ^(x) denotes the non-singlet Wilson coefficient with p qq (x) = 2/(1 — x) + — (1 + x). 
The singlet Wilson coefficients are of a similar structure. The splitting function P qq {x) related to the 



color factor d a j, c d ( I7.5|7.6p behaves like (1 



x 



,k>l 



ln(l — x) for large values of x, [391] . 



In Mellin-space the resummed Wilson coefficients have the structure 



C(N) = g (a s ) exp {ln(i%i(A) + g 2 (X) + a 8 g 3 (X) + 0(a 2 J(X))} 



14.13) 



with go{a s ) the normalization and A = a s /?oln(iV). In the flavor non-singlet case the functions gi, 2 (a s 
are given by 



gi{a s ) 

92(a s ) 



A, 



[A 



jeAi - Bi 
A) 



A)ln(l-A)] 

ln(l-A)-^[A + ln(l-A)] + ^ 

Po Po 



Aln(l-A) + -ln 2 (l-A) 
2 



with 



A x = 4CV, 1 2 



8C, 



/67 



[ Ts -c, c A -. Nl 



[393] fii = -3C F [394 



(14.14) 



(14.15) 



(14.16) 



The universal part of #3 oc a s ln(iV) was derived in [395] . 

The leading order resummation corrections beyond 0(a%) valid for large values of N in the non-singlet 
case have been calculated in [396] . Next-to-leading-log contributions were accounted for in [391(1397] . 

One may investigate the large- x structure of the 4-loop anomalous dimensions studying scheme- 
invariant evolution equations, cf. [15214155] . for the F 2 , F^-system, with <p a scalar particle coupling to 
the gluon as the photon couples to the quarks. The specific structure of the scheme- invariant kernels 
to the level of is assumed to hold at a*. Knowing the 3-loop Wilson coefficients, one may determine 
under this assumption the large x behaviour of the singlet splitting functions P^ (x) in the highest two 
powers in ln(l — x), cf. [398] . A resummation of the contributions to P qg , P gq , C 2t9 and Ci t9 has been 
performed in [392] . A systematic approach to obtain even higher order terms in case of the large-x 
resummation has been proposed in Ref. [399J. 
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15 Sum Rules and Integral Relations 



Deep-inelastic structure functions obey a series of sum rules for special moments or even integral re- 
lations, which are of interest for experimental tests. These relations are of different rigor. In most of 
the cases they receive radiative corrections, mass corrections, and in some cases even non-perturbative 
corrections a priori. Using current algebra techniques many sum rules have been investigated in [400J. 
A larger series of sum rules having been proposed for polarized scattering have been analyzed in [94J . 
Quark mass and QCD corrections have been given in [2BB] I400] . In the following we give a brief descrip- 
tion of the main sum rules and comment on the status of available QCD corrections. 

Adler sum rule |401j : 

This sum rule is rigorous and neither obtains QCD nor mass corrections. 

f 1 dr 

/ - [F?(x, Q 2 ) - F?(x, Q 2 )} = K(Nf), (15.1) 
Jo x 

with if (3) = 2 + 2sin 2 # c , (SU F (3)) and K(4) = 2, (SU F (A)), and 6 C the Cabibbo angle, cf [288M00] . 
Unpolarized Bjorken sum rule |402] : 

The sum rule refers to the charged current structure functions F± 

f dx (*, Q 2 ) - F?(x, Q 2 )] = K(Nf)A Fl (N f , Q 2 ), (15.2) 
Jo 

with if (3(4)) = 1 + sin 2 C , (1) and A F (N f , Q 2 ) = 1 + 0(a s ). The 3-loop QCD corrections to A Fl were 
given in [403| . 

Gross-Llewellyn Smith sum rule [404J: 

Likewise, the combination for the charged current structure functions F 3 yields 

f dx [F^(x, Q 2 ) + F» p (x, Q 2 )] = K(Nf)A Fa (N f , Q 2 ), (15.3) 
Jo 

with if (3(4)) = 6-2 sin 2 C , (6). The 3-loop QCD corrections to A Fs (N f , Q 2 ) were given in [305] and 
the 4- loop corrections in |406j . 
Polarized Bjorken sum rule [407J: 

The sum rule refers to the flavor non-singlet combination 

j'dx [gT(x,Q 2 )-gr(x,Q 2 )}=^ 

with gAy the neutron decay constants, QaI 9v ~ —1.26. The 3-loop QCD corrections to A gl (Nf, Q 2 ) 
were given in |4U5] and the 4-loop corrections in |4U8] . 
Gerasimov-Drell-Hearn sum rule [40911410] : 

This sum rule is given by the first moment of the polarized structure function g^' n (x,Q 2 ) in the form, 
cf. e.g. [338HHI1I1I2], 

M 2 f xo f --u 2 O 2 -> 

I P AQ 2 ) = 2^ J o dx 9 r(x, Q 2 ) = [ 2 ^n^ J _, OO (15 ' 5) 

at proton and neutron targets, with x = Q 2 / (2Mm n + m 2 +Q 2 ), fi Pt n the anomalous magnetic moment 
of the proton or nucleon (jl.ip . and Ti the first moment of the structure function gi at infinite space-like 
momentum transfer. The sum-rule has a very strong Q 2 -evolution for low values of the virtuality. In 
case of proton targets it changes sign between Q 2 = and Q 2 1 GeV 2 , [338]. 



gA 

gv 



A Sl (Nf,Q 2 ), 



(15.4) 
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Burkhardt-Cottingham sum rule [413J: 



1 

2 



dxg 2 (x,Q 2 ) = . (15.6) 



In the derivation of the sum rule it was assumed that the large //-behaviour of the associated amplitude 
A 2 , cf. [94], is governed by Regge-theory and it behaves as A 2 (Q 2 , v) ~ ^ _1_e , e < 0. However, this may 
be questioned, see |414j . It was argued in |415] that Regge cuts spoil the sum rule, which would vanish 
for large Q 2 , however. The sum-rule may be invalidated in case of a short-distance singularity |414j . see 
also |416j . The Burkhardt-Cottingham sum rule cannot be expressed in terms of expectation values of 
(axial-)vector current operators and not be derived using the light-cone expansion, see e.g. [94], although 
a formal analytic continuation for general values of the Mellin variable iV would suggest it. It is well- 
known, however, that specific moments may have a particular behaviour. The 0(a s ) corrections [400J , 
target mass corrections [93J and massive quark corrections |417] do not alter this relation. A brief review 
on this sum rule has been given in Ref. [414] . 
Efremov-Teryaev- Leader sum rule [418] : 

This sum rule refers to the valence (V) contributions of the polarized structure functions g\p related to 
the twist-3 operator matrix element d\ q ' v 



1 , e 2 



dx x 







9 r v (x } Q 2 ) + 2gJ^(x, Q 2 ) = = . (15.7) 



It receives quark mass corrections and is thus only valid in the limit m q — > 0, cf. [94J. It holds under 

certain conditions also in the presence of target mass corrections [93] . 

Ellis-Jaffe sum rule [419] : 

This sum rule is given by the integral 

f dxgl {n \x,Q 2 ) = C NS (Q 2 )A NS (Q 2 ) + C S (Q 2 )A S (Q 2 ). (15.8) 
Jo 

It was originally devised to describe the nucleon spin composition. However, the nucleon spin receives 
also gluonic and angular momentum contributions. The latter can in principle be measured using non- 
forward scattering [SO] ■ Here C NS ' S denote Wilson coefficients and A NS,S the corresponding OMEs in the 
flavor non-singlet and singlet case. The non-singlet current is conserved for massless quarks, while the 
singlet current is not conserved, which causes a non- vanishing singlet anomalous dimension at 0(a 2 ). 
The 3- and 4-loop QCD corrections to this sum rule have been calculated in [420(1421] . 
Gottfried sum rule [422] : 



f ^[F™(x,Q 2 )-F2 n (x,Q 2 )] = ~ [ dx[u v (x,Q 2 )-d v (x,Q 2 )-(d(x,Q 2 )-u(x,Q 2 )))} .(15.9) 
Jo x Jo 

This sum rule is non-rigorous since the distributions d(x, Q 2 ) and u(x, Q 2 ) are found to be different and 
are non-perturbative quantities. It is currently studied in lattice simulations, cf. Figure [TBI 
Integral Relations at Twist 2 and 3 [9311266] : 

The twist- 2 and twist-3 contributions to the structure functions gi,2{x,Q 2 ) obey the relations 



9 r\x,Q 2 ) = -g T r 2 (x,Q 2 )+ f 1 %Z =2 (v,Q 2 ) 

J x y 

9r 3 & Q 2 ) = Q 2 ) - 2 T %r\y, Q 2 ) 

L Jx v 



Q 2 



(15.10) 
(15.11) 
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The Wandzura-Wilczek relation (115. 10p [266J, originally derived in the massless case, remains valid in 
case of quark- and target mass corrections [93J, for the gluonic contributions [265], for non- forward scat- 
tering [301 ,423j , and diffractive scattering [30 211424] . It can be derived using the covariant parton model, 
as well as relations for other polarized structure functions containing twist-3 contributions, cf. [145,150j. 
In the limit of vanishing target masses g\ does not receive a twist-3 contribution. Eqs. (I15.10|ll5.1ip 
allow to disentangle the fermionic twist-2 and 3 contributions to g± and #2- 



16 Higher Twist Contributions 

Beyond the leading twist contributions, i.e. twist r = 2 in Quantum Chromodynamics, a quite different 
picture emerges for the structure functions both in the polarized and unpolarized case. With the twist-3 
contributions to polarized deep-inelastic scattering the problem of the operator mixing becomes more 
involved since the number of operators contributing at total spin N is growing with N. Both the 
anomalous dimensions and Wilson coefficients have been calculated during the last three decades at the 
level of the 0(a s ) corrections. The twist-3 anomalous dimensions have been derived in Refs. |425[ r426j . 
The evolution in the non-singlet case in the large N c limit has been studied in |427] and terms beyond 
this limit were derived in [428J. The twist-3 non-singlet and singlet anomalous dimensions have been 
calculated using light-ray operators in |429[I430] . The leading order Wilson coefficients were computed 
in |431] in the non-singlet and in |432] in the singlet case. The renormalization of gauge invariant 
operators contributing to the structure function g 2 (x,Q 2 ) has been studied in |433] . In Ref. |434] the 
gluonic contributions to the structure function g2(x, Q 2 ) were calculated. The complete non-singlet and 
singlet evolution of the twist-3 moments was derived in [435J. 

The data on g2(x,Q 2 ) |122[I135] have still large errors, which can be substantially improved at a 
facility like the EIC |137[ I138]. For first estimates on the behaviour of the twist-3 part of g2{x,Q 2 ) 
see e.g. in [436] . Yet the experimental errors on the twist-3 contribution of the structure function 
g%(x, Q 2 ) are large [295J. First lattice measurements of the lowest twist-3 moment have been performed 
in |43 7(1438] at larger pion masses and have to be refined for realistic values of in the future. 

In the unpolarized case the contributions of higher (dynamical) twist are of r = 4,6, .... They are 
largely suppressed in the limit of large virtualities Q 2 . However, the experimental data often exhibit a 
correlation between x and Q 2 due to similar values of S. Furthermore, in neutral current deep-inelastic 
scattering the largest statistics is located in the region of lower values of Q 2 . It is difficult to decide from 
which scale Q 2 onwards a data sample is widely free of higher twist contributions, for which sometimes 
a cut of W 2 = 12.5 — 15GeV is proposed [199,439j. This cut has been verified in the non-singlet 
analysis |199] . 

The local operators of higher twist can be constructed systematically near the light-cone. They are 
formed by more external quark and gluon fields than the twist-2 operators and potential contributions 
of lower twist operators in case mass scales are present. A twist-4 operator is given e.g. by 

: i){x)-i^d^ 2 ...d^ m i)(x)i){y)-i Vl d U2 ...d Un il){y) : (16.1) 

A systematic twist-decomposition has to be performed, cf. [140J . One forms operator matrix elements 
with these operators between nucleon states. A 'partonic' interpretation assumes, that all external 
lines can be factorized. The corresponding contributions to deep-inelastic structure functions are then 
similar as in the case of twist-2, with the generalization that both the Wilson coefficients and OMEs 
depend on n — 1 dimensionless parameters Y17=i Xi = x i wnere n denotes the number of external fields. 
These parameters, except x, cannot be measured in the deep-inelastic process. 
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Figure 22: Comparison of the higher twist coefficient Cht{x) = CJ =4 (x) in the large x region for the proton data 
as function of x in a NLO (dotted lines), NNLO (dashed lines), N3LO analysis (dash-dotted lines) for the non-singlet 
QCD Wilson coefficient (full lines). Some bin centers are slightly shifted for better visibility; from |454j, ©(2010) 
by Elsevier Science. 



Early theoretical investigations of the structure of higher twist operators and their anomalous di- 
mensions in D = 6 $ 3 -theory [440J and QCD [186,425 144114447] revealed the principal structure of these 
contributions. The lowest order anomalous dimensions have been calculated in |186[l444ll447] and the 
Wilson coefficients, in different operator bases, in |442ll443ll445] . More recently also gluonic operators 
were considered |446] . A systematic study of higher twist light-cone distribution amplitudes was given 
in |447j . The renormalization of these operators has been worked out systematically in |448] . The 
evolution of the lowest twist-4 moments at leading order has been studied in |449] recently. 

Due to the fact that the OMEs of higher twist operators cannot be measured experimentally one 
may envisage their determination with lattice simulations. At present, investigations of this kind are 
still at the beginning. In Mellin-space the basis of higher twist operators grows considerably with the 
value of the Mellin variable iV and an according number of OMEs has to be measured on the lattice, 
which further complicates quantifications of higher twist effects ab initio. 

Estimates of higher twist effects have been obtained studying renormalon corrections to sum-rules 
and deep-inelastic structure functions, see |450] and for reviews Refs. [4511 P452J . 

Because of the complications mentioned the extraction of higher twist contributions to deep-inelastic 
structure functions is usually being performed applying a suitable Ansatz for these contributions: 



F^Q 2 ) = Fr 2 (x,Q 2 



16.21 



Q 2 Q 2 

Sometimes one considers the higher twist terms in (116. 2p in additive form. In many cases these con- 
tributions are fitted together with the PDFs, cf. [453, 454J . Since the logarithmic scaling violations of 
the higher twist coefficients are in general different from those of the lowest twist it was suggested to 
proceed in a different way, cf. |199[l454j . In the flavor non-singlet case the lowest twist contributions to 
the structure functions may be determined in the region which is free of higher twist terms. One then 
extrapolates the twist-2 distributions into the region, where higher twist contributions are present and 
measures the difference. The size of higher twist contributions thus obtained stabilizes including higher 
order QCD corrections to the leading twist terms. The inclusion of the leading twist 3-loop Wilson 
coefficients yields the dominant contribution and effectively describes the N 3 LO terms, while the N 4 LO 
term is approximate only as it refers to the soft exponentiations |396j . In Figure [221 the corresponding 
results are shown for a world analysis in the valence region for proton and deuteron data. 
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Figure 23: Left: comparison between the higher twist values corresponding to the fits of inclusive DIS LSS06 [455J 
and combined inclusive and SI DIS data set LSS10 [342| ©(2010) by the American Physical Society. Right: The 
additive higher twist coefficients Cp(x),Cd(x) and C n (x) [329], ©(2010) by Elsevier Science. 

A similar analysis is not possible in the polarized case at present, because of the smaller kinematic 
range in Q 2 and the larger errors of the data for the structure function g\{x, Q 2 ). Here the higher twist 
contributions are fitted as additive terms in an NLO analysis for individual bins, cf. e.g. [32911342, 455]. 
In Refs. [34211455] smaller errors have been obtained with an excess at lower values of x in case of 
neutron targets and proton targets, see Figure [231 In [329] higher twist contributions which are widely 
compatible with zero have been obtained both for the proton and neutron targets, with a 1-1.5 a effect 
in two bins for the proton data. To reveal twist-4 contributions to the structure function gi(x,Q 2 ) 
more precise data from future experiments are needed. 



17 Nuclear Parton Distribution Functions 

The parton densities for free nucleons and nuclei are different. This has been first demonstrated in detail 
by the EMC-experiment [456]. One well-known effect is due to the Fermi-motion of nucleons within 
nuclei [457] . Beyond these contributions other deviations have been observed, in the region of lower 
values of x. Various phenomenological explanations have been proposed, like binding models of nucleons, 
admixture models, change of nucleon mass, pion enhancement, multi-quark clusters, dynamical rescaling 
(change of confinement size), shadowing, and others; for reviews see [458] . Nuclear effects also modify the 
Drell-Yan process, W^/Z-production, highpx hadron production, J/ip- and direct photon production. 

Recent global data analyses to extract the nuclear PDFs have been carried out by different groups 
[459 d464j extracting nuclear correction factors Rf(xN,Q 2 ) for the parton densities fi(x,Q 2 ) in free 
nucleons. In general the Bjorken variable x^ for nuclear PDFs has support x^ G [0, A]. However, 
measured distributions for values of x^ > 1 do rapidly fall. Therefore, one often parameterizes the 
nuclear PDFs only for the range x?j G [0, 1]. A recent parameterization was given in Ref. [464] . The 
correction factors Rf(xN, Q 2 ) are illustrated in Figure 1241 for a series of nuclei for the quark distributions 
u v ,u, s and the gluon distribution g. The d v distribution is assigned the same correction factor as the 
u v distribution. 

The nuclear corrections at very small and very large x^ become stronger with the rising value of 
A. At large xn Fermi- motion effects dominate. Anti-shadowing effects are observed in the region 
xn G [0.01, 0.1] and the correction factors become smaller than one for even smaller values of xn- The 
comparison of different fits shows that there are still quite some uncertainties. On the theoretical side 
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a complete understanding of the factors Rf is not yet obtained, despite a large number of proposed 
models. 




Figure 24: The NLO nuclear modification factors Rf for different parton densities at Q 2 = lOGeV 2 for Be, Fe, Au 
and Pb targets as a function of x^. The inner (outer) error bands correspond to Ax 2 = 1(30), respectively. Red 
lines: |464j : dash-dotted lines |463| : dashed lines |460| : from |464j ©(2012) by the American Physical Society. 

Deep- inelastic data analyses include deuteron data to obtain an improved representation of the down- 
quark distributions due to scattering off the neutron. The deuteron data require nuclear corrections 
[465J, which mainly concern Fermi-motion and the off-shell effect, cf. [117] . For the off-shell corrections 
the Bonn [466] and Paris |467] potentials are used, including the induced differences into the systematic 
errors. One assumes that the nuclear model suggested in [461] can be applied to the case of light nuclei, 
like deuterium, which has recently been confirmed for 3 He and 9 Be targets [462] . Di-muon production 
data, which are important to measure the strange-quark distributions, are currently measured at nuclear 
targets, for which the corrections [468] were applied in [117] . Drell-Yan data of nuclear targets which 
are used in present PDF- fits are measured in a region with lower nuclear effects, cf. [469J. 

18 Conclusions 

During the last 40 years the understanding of the sub-structure of nucleons took a vast development. 
The theory of strong interactions, QCD, can now be probed at the per-cent level in an increasing 
number of reactions. High-luminosity measurements performed at SLAC, CERN, DESY, FERMILAB, 
and JLAB more and more revealed the dynamics of quarks and gluons at short distances. With the 
advent of the LHC this process will continue in an intimate partnership between precision measurements 
and precision calculations. During the coming years the data analysis of the HERA-experiments will 
be finalized. Here a remaining issue is the theoretical understanding of the ep — > jet cross sections at 
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NNLO. Jet physics, the production of W ± and Z-bosons, the Drell-Yan cross section, and inclusive 
tt- as well as single-top production will be important scattering cross sections measured at the LHC 
to limit the present experimental errors on a a (M§) and to further refine the knowledge on the parton 
distributions. The precision of the data requests NNLO QCD corrections. 

In the more distant future high luminosity facilities like EIC may allow to answer many open 
questions related to lower energies, like the partonic structure of polarized nucleons, twist-3 effects, 
and higher twist effects in general, and unravel transverse degrees of freedom, including associated spin 
effects in detail. This machine will also have the potential to measure a s (M|) very precisely. Proposed 
facilities like LHeC, on the other hand, allow to further progress into the small x and high Q 2 domain 
by ep scattering. Many proposed theoretical concepts can be tested in this way at high precision, giving 
a boost to even more refined theoretical calculations than carried out at present. Probing ever shorter 
distances may finally answer the more fundamental question for a possible substructures of quarks, 
which will request precision, both at the experimental and theoretical side. 
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